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(57) Abstract: The present invention identifies 
dififerent activities mediated by eosinophil PGD2 
receptors and features methods measuring the 
ability of a compound to modulate such activities. 
Activities mediated by eosinophil PGD2 receptors 
include those associated with CRHT2 and those 
associated with the DP receptor. Activities 
identified herein as associated with eosinophil 
CRHT2 include a change in cell morphology, 
degranulation, and a specific chemokinetic effect: 
Activities identified herein as associated with 
the eosinophil DP receptor include resistance to 
apoptosis. 
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TITLE OF THE INVENTION 

EOSINOPHIL PROSTAGLANDIN D2 RECEPTOR ASSAYS 

CROSS-REFERENCE TO RELATED APPLICATIONS 
5 The present application claims priority to provisional application U.S. 

Serial No, 60/306,357, filed July 18, 2001, hereby incorporated by reference herein. 

BACKGROUND OF THE INVENTION 

The Background of the Invention and references cited in the present 
10 application are not admitted to be prior art to the claimed invention. 

Prostanglandin D2 (PGD2) is a cyclooxygenase metabolite of 

arachidonic acid. (Narumiya, etaL, Physiological Reviews 79:1193-1226, 1999.) 
PGD2 has been implicated in playing a role in different physiological events such as 

sleep and allergic responses. (Boie, et al, The Journal of Biological Chemistry, 
15 270:18910-18916, 1995, Narumiya, et al. Physiological Reviews 79:1193-1226, 
1999, Matsuoka, et al. Science 257:2013-2017, 2000.) 

Mast cells and TH2 cells are important inmiune cells involved in 
allergic responses. PGD2 is released from mast and TH2 cells in response to an 

immunological challenge. (Roberts, et aL, N. Engl J. Med. 505:1400, 1980, Lewis, et 
20 a/., /. Immunol 129:1621, 1982, Tanaka, et al, J Immunol 164:2211, 2000.) 

Receptors for PGD2 include the "DP" receptor, the chemoattractant 

receptor-homologous molecule expressed on TH2 cells ("CRTH2"), and the "FP" 
receptor. These receptors are G-protein coupled receptors activated by PGD2. PGD2 

is a non-selective agonist at the FP receptor. (Abramovitz, et al, Biochimica et 
25 Biophysica Acta 1483:2S5'293, 2000.) 

Abramovitz, et al^ U.S. Patent No. 5,958,723 and Boie, et al.. Journal 
of Biological Chemistry 270:18910-18916, 1995, describe the cloning and 
characterization of the human DP receptor. These references also indicate that PGD2 

activates the DP receptor. 
30 Abe, et al. Gene 227:11-11, 1999, Nagata, et al, FEBS Letters 

459:195-199, 1999, and Nagata, et al.. The Journal of Immunology 7(52:1278-1286, 
1999, describe CRTH2 and its expression on different cells including human T-helper 
cells, basophils, and eosinophils. Hirai, et al, J, Exp. Med. 795:255-261, 2001, 
indicates that CRTH2 is a receptor for PGD2. 
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SUMMARY OF THE INVENTION 

The present invention identifies different activities mediated by 
eosinophil PGD2 receptors and features methods measuring the ability of a compound 
to modulate such activities. Activities mediated by eosinophil PGD2 receptors 
5 include those associated with CRHT2 and those associated with the DP receptor. 
Activities identified herein as associated v^th eosinophil CRHT2 include a change in 
cell morphology, degranulation, and a specific chemokinetic effect. Activities 
identified herein as associated with the eosinophil DP receptor include resistance to 
apoptosis. 

10 Measuring the ability of a conipound to modulate a PGD2 receptor 

activity can be performed quantitatively or qualitatively. Compounds modulating 
PGD2 receptor activity include agonists, antagonists and allosteric modulators. 

Thus, a first aspect of the present invention features a method that 
measures the effect of a test compound on either apoptosis or degranulation as a 
15 measure of the ability of the compound to modulate a PGD2 receptor activity. The 

method employs eosinophil cells. 

Another aspect of the present invention describes a method of assaying 
the ability of a test compound to modulate CRTH2 activity using a compound 
identified as binding to GRTH2. The method comprises the steps of; (a) identifying a 

20 compound that binds to human CRTH2; (b) providing the compound to an eosinophil; 
and (c) measuring eosinophil morphology, chemokinesis under conditions 
distinguishing chemokinesis from chemotactic ability, or degranulation. 

Another aspect of the present invention describes a method of assaying 
the ability of a test compound to modulate CRTH2 activity involving the use of a 

25 CRTH2 agonist. The method comprises the steps of: (a) providing the test compound 
and an CRTH2 agonist to an eosinophil, and (b) measuring eithei: eosinophil 
niorphology, chemokinesis under conditions distinguishing chemokinesis from 
chemotactic ability, or degranulation. 

Other features and advantages of the present invention are apparent 

30 from the additional descriptions provided herein including the different examples. 
The provided examples illustrate different components and nietiiodology useful in 
practicing the present invention. T!he examples do not limit the claimed invention. 
Based on the present disclosure the skilled artisan can identify and employ other 
components and methodology useful for practicing the present invention. 

35 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates PGD2 receptor expression on human eosinophils. 

RT-PCR was performed on total RNA isolated from purified human eosinophils using 
CRTH2 or DP-specific primers. The RT-PCR product was revealed, by Southern 
5 blot, using a CRTH2-specific (top panel) or DP-specific Qower panel) radioactive 
probe. RNA from HEK cells expressing recombinant CRTH2 or DP receptor was 
used as a positive control in lane 1. RT-PCR using 18S ribosomal RNA-specific 
primers was conducted in pairallel to ensure that equivalent amounts of RNA were 
used between each donor (data not shown). The bands seen are derived from mRNA 
10 and not genomic DNA since no signal is detected in absence of reverse transcriptase. 
Results from two out of four donors tested are shown. 

Figure 2 illustrates a rapid change in eosinophil morphology induced 
by PGD2. Purified human eosinophils were incubated for 15 minutes with various 

agents in a 24-well dish. Cells were then magnified 200-times using an inverted 
15 microscope, a, vehicle-treated eosinophils, b, eosinophils treated with 10 nM PGD2. 
c, 1 |iM BW245C (a DP-selective agonist) d, 10 nM 13,14-dihydro-15-keto-PGD2 
(DK- PGD2). e, 100 nM platelet-activating factor; PAF. f, 1 ng/ml of interleukin-5; 

IL-5. A representative experiment from 20 donors tested is shown. 

Figure 3 illustrates the effect of PGD2 on eosinophil chemokinesis. 

20 Purified human eosinophils were treated for 5 minutes with various agents prior to 
being placed in the upper chamber of a chemotactic unit. No chemoattractant was 
added to the lower chamber in order to simply measure chemokinesis. After two 
hours, the number of cells that transmigrated to the lower chamber was evaluated with 
an hematocytometer. Chemokinesis efficiency is expressed as the number of 

25 transmigrating cells with the agent divided by the number of transmigrating cells with 
vehicle only (fold-increase chemokinesis over biackground). Lane 1, vehicle treated 
eosinophils. Lane 2, eosinophils were treated with 100 nM PGD2, lane 3 with 1 |JiM 
BW245C, lane 4 with 100 nM DK- PGD2, lane 5 with 100 nM of platelet activating 
factor, lane 6 with 1 ng/ml of interleukin-5 and lanes 7-8-9 with 1 \jM of the indicated 

30 compounds. For each experiment, each condition was tested in two independent 
wells. The mean response is indicated by a dash. The effect of PGD2 at 100 nM is 

significant with a probability of < 0.001 in repeated measures ANOVA followed by 
paired t-tests. 

Figure 4 illustrates the ability of PGD2 to trigger eosinophil 
35 degranulation. Purified human eosinophils were treated for 1 hour with various 
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agents. The amount of ECP released in the media was then determined by 
radioimmunoassay. Lane 1, vehicle treated cells. Lane 2, eosinophils were treated 
with 100 nM PGD2, lane 3 with 1 pM BW245C, lane 4 with 100 nM DK-PGD2, lane 
5 with 100 nM of platelet activating factor, lane 6 with 1 ng/ml of interleukin-5 and , 

5 lanes 7-8-9 with 1 |xM of the indicated compounds. The values correspond to the 
amount of ECP detected under the various conditions minus the value obtained with 
the vehicle treated cells (amount of ECP over background). For each experinient, 
each condition was tested in duplicate. The mean response is indicated by a dash. 
The effect of PGD2 at 100 nM was tested on 11 donors (p<0.0003 in t-test) while the 

10 effect of DK- PGD2 at 100 nM was tested on 8 donors (p<0.01). The value presented 

for PAF is the mean of six independent experiments (p<0.dl). 

Figure 5 illustrates the ability of PGD2 to increase the survival of 

eosinophils in culture. Purified human eosinophils were maintained in culture in the 
presence of various agents for 36 hours. The cells were then harvested and the extent 

15 of apoptosis was evaluated by flow cytometry (Annexin V/propidium iodide staining). 
Cells that have not reached the stage of late apoptosis (thus not positive for both 
annexin V and propidium iodide staining) were considered to be alive. Lane 1, 
vehicle treated cells. Lane 2-9, eosinophils treated with 1 nM of the indicated 
compounds except lane 6 where interleukin-5 was used at 1 ng/ml. The values 

20 correspond to the percentage of non-late apoptotic eosinophils in the treated 

population minus the percentage of non-late apoptotic eosinophils in the vehicle 
treated population. The mean response is indicated by a dash. The effect of PGD2 at 
1 (jM was tested on 7 donors (p<0.1 in t-test) while the effect of BW245C at 1 \M 
was tested on 9 donors (p<0.0005). 

25 

DETAILEP DESCRIPTION OF THE INVENTION 

Identifying different effects mediated by eosinophil PGD2 receptor 

activation provides for indicators that may be measured to evaluate the ability of a 

compound to modulate eosinophil PGD2 receptor activity and provides information 
30 concerning the physiological effects of PGD2 receptor activation. Information 

concerning the physiological effects of PGD2 receptor activation can be used to help 

evaluate the importance of inhibiting a PGD2 receptor activity. 

Compounds modulating eosinophil PGD2 receptor activity have a 

variety of different uses including utility as a tool to further study PGD2 receptor 
35 activity and as an agent to achieve a beneficial effect in a patient. Modulating PGD2 
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receptor activity includes evoking a response at the receptor and altering a response 
evoked by a PGD2 receptor agonist or antagonist. 

Beneficial effects of modulating PGD2 receptor activity include 

achieving one or more of the following in a patient: the treatment or prevention of an 
5 inflammatory disease such as asthma, treatment or prevention of allergic rhinitis or 
arthritis; and the treatment or prevention of a sleep disorder. A patient is a mammal, 
preferably a human. Reference to patient does not necessarily indicate the presence of 
a disease or disorder. The term patient includes subjects treated prophylactically and 
subjects afflicted with a disease or disorder. 
10 Selective agonists or antagonists that mimic or block PGD2 actions at 

the DP receptor, CRTH2 and/or FP receptor may have utility in the treatment of 
disease states or diseases including but not limited to allergic rhinitis and other 
allergic conditions in which mast cells, eosinophils, TH2 cells and other inraiune cells 
express the DP receptor, CRTH2, and/or FP receptor, or produce PGD2. Additional 

15 examples of therapeutic applications include one or more of the following: sleep 
disorders; glaucoma; osteoporosis; modulators may be useful as cytoprotective, 
analgesic or anti-inflammatory agents; modulators inhibiting platelet aggregation may 
be useful for treating vascular disease, prevention of post-injury blood clotting, 
rejection in organ transplant and by-pass surgery, congestive heart failure, pulmonary 

20 hypotension and Raynaud's disease. 

Eosinophil PGD^ Receptor and Inflammation 

Eosinophils were found to express the DP receptor and CRTH2. The 
different effects mediated by PGD2 at these receptors appear to assist the 
25 inflanamation response. Pharmacological blockade of PGD2-mediated events at both 
the DP receptor and CRTH2 may reduce damage caused by eosinophils at an 

inflammation site. 

In allergic situations, PGD2 is released by mast cells and may facilitate 

entry into the inflammation site through DP-mediated vasodilation/extravasation of 
30 eosinophils as well as other circulating leukocytes. (Mantovani, et ah. Lancet 

545:1499, 1994). The entry of eosinophils into the allergic site would be stimulated 
by the pro-chemokinetic activity of PGD2 through CRTH2. 

The anti-apoptotic and degranulation effects of PGD2 on eosinophils 
appear to be playing a factor in inflammation. The survival of resident eosinophils 
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would be prolonged by the anti-apoptotic effect of PGD2 acting through the DP 
receptor. 

Eosinophil degranulation triggered by PGD2 activation through 
CRTH2 causes the release of granule-derived proteins. The effects of granule 
5 proteins include cytotoxicity at the bronchial epithelium, an increase in nonspecific 
bronchial hyperreactivity and impaired ciliary function. 

PGD? Receptor Assays 

Different types of assays formats can be employed making usie of the 
10 activities identified herein as associated with the eosinophil DP receptor or the 
eosinophil CRHT2. Examples of such formats include: 

(1) Measuring the ability of a compound to affect apoptosis or 

degranulation; 

(2) Identifying a compound that binds to the human eosinophil 

15 CRHT2 and then testing the ability of the compound to affect eosinophil morphology, 
chemokinesis under conditions distinguishing chemokinesis from chemotactic ability, 
or degranulation; and 

(3) Screening for a CRHT2 antagonist using a CRHT2 agonist and 
measuring the ability of a test compound to modulate changes in eosinophil 

20 morphology, chemokinesis under conditions distinguishing chemokinesis from 
chemotactic ability, or degranulation produced by the agonist. 

Measuring the effect of a compound oh apoptosis or degranulation 
provides an overall measure of the effect of the compound on DP receptor or CRTEI2 
activity. Measuring apoptosis or degranulation also provides a direct measure on 

25 activities that it would be desirable to inhibit. 

Li an embodiment of the present invention, a binding assay is 
employed to select for compound binding to a prostaglandin D2 receptor prior to ah 

apoptosis or degranulation assay. Assays measuring the ability of a compound to bind 
to a DP receptor or CRTH2, employ a DP receptor or CRTH2 polypeptide comprising 
30 a PGD2 binding site. DP receptor and CRTH2 polypeptides include full-length 
human receptors and functional derivatives thereof, fragments containing a PGD2 

binding site, and chimeric polypeptides comprising such fragments. A chimeric 
polypeptide comprising a fragment that binds PGD2 also contains one or more 

polypeptide regions not found in a human DP receptor or CRTH2. 
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Preferably, assays measuring PGD2 binding employ full length human 
DP receptor or CRTH2. The human DP receptor is described by Abramovitz, et al. 
V S. Patent No, 5,958,723. Human CRTH2 is described by Nagata, et al, The 
Journal of Immunology 762:1278-1286, 1999, and Gen-Bank Accession No. 
5 AB00535. 

PGD2 receptor amino acid sequences involved in PGD2 binding can 
be identified using labeled PGD2 and different PGD2 receptor fragments. Different 
strategies can be employed to select fragments to be tested to narrow down the 
binding region. Examples of such strategies include testing consecutive fragments 
10 about 15 amino acids in length starting at the N-terminus, and testing longer length 
fragments. If longer length fragments are tested, a fragment binding PGD2 can be 
subdivided or mutated to further locate the PGD2 binding region. Fragments used for 

binding studies can be generated using recombinant nucleic acid techniques. 

Binding assays can be performed using recombinantly produced PGD2 

15 receptor polypeptides present in different environments. Such environments include, 
for example, cell extracts and purified cell extracts containing a PGD2 receptor 

polypeptide expressed from recombiniant nucleic acid or naturally occurring nucleic 
acid and also include, for example, the use of a purified PGD2 receptor polypeptide 

produced by recombinant means or from naturally occurring nucleic acid which is 

20 introduced into a different environment. 

The ability of a coinpound to antagonize PGP2 receptor activity can be 
evaluated using a PGD2 agonist able to produce receptor activity and then measuring 
the ability of one or more test compounds to alter such activity. Agonists that can be 
employed include those able to stimulate both DP receptor activity and CRHT2 

25 activity and those selective for DP receptor activity or CRHT2 activity. Examples of 
different types of agonists are PGD2 which acts at both die DP receptor and CRHT2; 
13-14-dihydro-15-keto-PGD2 which is specific for CRTH2; and BW245C which is 

specific for the DP receptor. 

The effectiveness of an antagonist to alter PGD2 receptor activity can 
30 be evaluated by comparing PGD2 receptor activity in the presence of the agonist with 
such activity in the presence of the agonist and antagonist. Different types of assay 
formats can be employed. For example, a control experiment involving an agonist 
and a test experiment involving the agonist and a test compound can be performed at 
the same or at different times. 
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Techniques for measuring apoptosis, morphology, chemokinesis under 
conditions distinguishing cheniokinesis from chemotactic ability, and degranulation 
are well known in the art. Changes in morphology can be measured visually with the 
aid of a microscope, such as by scoring cells with irregular shapes. Techniques for 
5 measuring morphology include those described in the Examples provided below. 

Apoptosis is a type of cell death that is programmed by the cell. 
Techniques for measuring apoptosis include those described in the Examples provided 
below. 

Chemokinesis is an increase in cell mobility that is brought about by a 
10 reagent in the absence of chemical gradient. Techniques for measuring chemokinesis 
include those described in the Examples provided below. 

Degranulation results in the release of granule-derived proteins, such 
as the major basic protein, the eosinophil cationic protein, eosinophil-derived 
neurotoxin, and eosinophil peroxidase. Techniques for measuring degranulation 
15 include those described in the E?camples provided below. 

Dosing For Therapeutic Applications 

Guidelines for pharmaceutical administration in general are provided 
in, for example, Remington 's Pharmaceutical Sciences 1 Edition^ Ed. Gennaro, 
20 Mack Publishing, 1990, and Modem Pharmaceutics 2^^ Edition, Eds. Banker and 
Rhodes, Marcel Dekker, Inc., 1990, both of which iare hereby incorporated by 
reference herein. 

PGD2 receptor active compounds having appiropriate functional groups 
can be prepaii^ as acidic or base salts. Pharmaceutically acceptable salts (in the form 

25 of water- or oil-soluble or dispersible products) include conventional nonrtoxic salts 
or the quaternary anmionium salts that are fonned, from inorganic or organic 
acids or bases. Examples of such salts include acid addition salts such as acetate, 
adipate, alginate, aspartate, benzoate, benzenesulfonate, bisulfate, butyrate, citrate, 
camphorate, camphorsuifonate, cyclopentanepropionate, digluconate, dodecylsulfate, 

30 ethariesulfonate, fumarate, glucoheptanoate, glycerophosphate, hemisulfate, 
heptanoate, hexanoate, hydrochloride, hydrobrbmide, hydroiodide, 2- 
hydroxyethanesulfonate, lactate, maleate, methanesulfonate, 2-naphthalenesulfonate, 
nicotinate, oxalate, pamoate, pectinate, persulfate, 3-phenylpropionate, picrate, 
pivalate, propionate, succinate, tartrate, thiocyanate, tosylate, and undecanoate; and 

35 base salts such as ammonium salts, alkali metal salts such as sodium and potassium 
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salts, alkaline earth metal salts such as calcium and magnesium salts, salts with 

organic bases such as dicyclohexylamine salts, N-methyl-D-glucamine, and salts with 

amino acids such as arginine and lysine. 

PGD2 receptor active compounds can be administered using different 

5 routes including oral, nasal, by injection, and transmucosally. Active ingredients to be 
administered orally as a suspension can be prepared according to techniques well 
known in the art of pharmaceutical formulation and may contain microcrystalline 
cellulose for imparting bulk, alginic acid or sodium alginate as a suspending agent, 
methylcellulose as a viscosity enhancer, and sweeteners/flavoring agents. As 

10 immediate release tablets, these compositions may contain naicrdcrystalline cellulose, 
dicalcium phosphiate, starch, magnesium stearate and lactose and/or other excipients, 
binders, extenders, disintegrants, diluents and lubricants. 

When administered by nasal aerosol or inhalation, compositions can be 
prepared according to techniques well known in the art of pharmaceutical 

15 formulation. Such techniques can involve preparing solutions in saline, employing 
benzyl alcohol or other suitable preservatives, absorption promoters to enhance 
bioavailability, fluorocarbons, or other sdlubiliang or dispersing agents. 

Routes of administration include intravenous (both bolus and 
infusion), intraperitoneal, subcutaneous, topical with or without occlusion, and 

20 intramuscular. Injectable solutions or suspensions known in the art include suitable 
non-toxic, parenterally-acceptable diluents or solvents, such as mannitol, 1,3- 
butanediol, water. Ringer's solution and isotonic sodium chloride solution. Dispersing 
or wetting and suspending agents, include sterile, bland, fixed oils, such as synthetic 
mono- or diglycerides; and fatty acids, such as oleic acid. 

25 Rectal administration in the form of suppositories, include the use of a 

suitable non-initating excipient, such as cocoa butter, synthetic glyceride esters or 
polyethylene glycols. These excipients are solid at ordinary temperatures, but 
liquidify and/or dissolve in the rectal cavity to release the drug. 

Suitable dosing regimens for therapeutic applications can be obtained 

30 taking into account factors well known in the art including age, weight, sex and 

medical condition of the patient; the severity of the condition to be treated; the route 
of administration; the renal and hepatic function of the patient; and the particular 
compound employed. 

Optimal precision in achieving concentrations of drug within the range 

35 that yields efficacy without toxicity requires a regimen based on the kinetics of the 
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drug's availability to target sites. This involves a consideration of the distribution, 
equilibrium, and elimination of a drug. The daily dose for a patient is expected to be 
between 0.01 and 1,000 mg per adult patient per day. 

5 Examples 

Examples are provided below to further illustrate different features of 
the present invention. The examples also illustrate useful methodology for practicing 
the invention. These examples do not limit the claimed invention. 

10 Example 1: Material and Meithods 

This example illustrates different reagents and techniques. 

Reagents 

PGD2, fluprostenol and PGE2 were obtained from Biomol Research 
15 Laboratories, (Plymouth Meeting, PA). BW245C, 13,14-dihydro-15-keto-PGD2 
(DK- PGD2) and latanoprost (free acid) were from Cayman Chemical (Ann Arbor, 
MI). Platelet activating factor (PAF) was from Sigma (St-Louis, MO). Recombinant 
human interleukin-5 was produced using a baculovirus system and purified by FPLC. 
(Brown, et al. Protein Expr. Purif. (5:63, 1995.) 

20 

Eosinophil Purification 

Circulating eosinophils were isolated from heparinized venous blood 
from normal volunteers. Erythrocytes were removed by addition of Dextran to a final 
concentration of 0.9% (Dextran TSOO from Pharmacia prepared as a 6% stock in 0.9% 
25 saline solution). After a 45 minute incubation at room temperature, the leukocytes in 
the plasma fraction were collected by centriftigation (4 °C, 300xg, 10 minutes), and 
resuspended in Hank's balanced salt solution (BBSS without calcium and 
magnesium). 

A density step gradient was generated by placing 20 ml of Ficoll- 
30 Paque™ (Pharmacia) under 30 ml of resuspended cells. The gradient was centriftiged 
(4 °C, 400xg, 30 minutes) and the pellet containing die granulocytes was resuspended 
in 10 ml of water for 15 seconds to lyse any residual erytfirocytes. The hypotonic 
lysis was stopped by the addition of 40 ml of HBSS. 

The cells were then centriftiged (4 °C, 3pOxg, 10 minutes), washed 
35 once with 50 ml of HBSS and resuspended in Dulbecco phosphate buffer saline (PBS 
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without calcium and magnesium from GIBCO-BRL) at a concentration of 1 x lO' 
cells per ml. An equal volume of CD 16 magnetic beads ^Iteny Biotec) was added 
and incubated at 4 for 30 minutes. At the end of the incubation, the volume was 
brought to 1 ml with PBS (without Ca*^ and Mg*^) and applied to a CS separation 
5 column placed in the magnetic field of a MACS separator (Milteny Biotec). 

The CD16+ neutrophils were retained in the column while a >95% 
pure fraction of CD16- eosinophils eluted from the column. The purity of the 
eosinophil fraction was evaluated by flow cytometry (CELL-DYN 3700 System) 
based on size, complexity, granularity and Ibbularity. 

10 

RT-PCR and Southern Blot 

Total RNA was obtained from isolated eosinophils (>95% pure) by 

using a total RNA isolation kit (Rneasy kit, Qiagen) and treated with DNAse (Gibco- 
BRL) prior to reverse transcription (Gene Amp kit, Perkin Elmer). Amplification of 

15 DP receptor by PCR (Advantage GC kit, Clontech) used the following primers: DP 
sense, 5'-ACAACTCGTTGTGCCAAGCC (SEQ, ID. NO. 1); DP antisense, 5'- 
GCATCGCATAGAGGTTGCGC (SEQ. ID. NO. 2); CRTH2 sense, 5'- 
CTACAATGTGCTGCTCCTGAAC (SEQ. ID. NO. 3); CRTH2 antisense, 5'- 
CAGGTGAGCACGTAGAGCAC (SEQ. ID. NO. 4). The PCR reaction (50 fil) 

20 included a denaturation step (94 "C, 1 minute) and 35 cycles of PCR (94 "C, 30 
seconds; 55 °C, 30 seconds; 68 °C, 1 minutes). 

PCR reactions were electrophoresed in agarose gels and transferred to 
nylon N+Hybond membrane (Amersham). The blot was hybridized with a ^^P-labeled 
DNA fragment encoding the full-length hCRTH2 or hDP receptor in ExprcssHyb 

25 solution (Clontech) overnight at 68 °C. The blot was washed twice iii 2X SSC (at 
65 "C) and twice in 0.2X SSC (at 65 °C) for 30 minutes each. Results were revealed 
by autoradiography. 

In Situ Hybridization 

30 Freshly isolated eosinophils (2x10^) were layered onto a poly-D- 

lysine coated glass slide by centrifugation (Cytospin). Cells were then fixed in 4% 
paraformaldehyde solution prepared in PBS (pH of 7.4) for 20 minutes at room 
temporature. The slides were then processed as follows: 2 minutes in 3X PBS, 2 
times 2 minutes in IX PBS and incubations of 5 minutes in 50%, 70%, 95% and 

35 100% aqueous ethanol solutions. Slides were air dried and stored at -80 °C. 
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The slides were thawed to room temperature and washed for 5 minutes 
in diethylpyrocarbonate (DEPC)-treated water and twice in PBS. The sections were 
treated with 1.0 )ig/ral proteinase K in 100 mM Tris, pH 8.0, 50 mM EDTA for 10 
minutes at 37 °C and washed for 5 minutes in DEPC-treated water. The slides werei 
5 then washed in 0.1 M triethanolamine, pH 8.0 (TEA) for 5 minutes and washed again 
for 10 minutes in TEA with 0.25% acetic anhydride. Finally, the sections were 
washed twice for 5 minutes in 2x SSC. 

A 398 bp fragment representing the 5' terminal end of the human DP 
receptor cDNA was amplified by PGR and subcloned into the PGR n dual promoter 

10 vector (Invitrogen). The plasmid was linearized using either Xho I or Spe I and 

digoxigenin-labeled (DIG) riboprobes were synthesized using the DIG-RNA labeling 
kit from Boehringer Mannheim. The riboprobes were diluted in 75% hybridization 
buffer (75% formamide, 3x SSG, Ix Denhardt*s, 0.2 mg/ml yeast tRNA, 50 mM 
sodium phosphate, 10% dextran sulfate) and layered onto the cytospin slides. The 

15 slides were covered with parafilm and left to hybridize for 16 hours at 55 ""C in a 
humidified (75% formamide) chamber. The parafilm was then removed by soaking 
the slides in 2x SSG for 30 minutes. The sections were then treated with RNase A (40 
|ig/ml in 10 mM tris, pH 8.0, 500 mM NaGl) for 45 minutes at 37 X. The slides 
were washed in 2x SSC, Ix SSG, 0.5x SSG (for 10 minutes each at room temperature) 

20 and b.lx SSC (45 minutes at 60°G). 

Colorimetric detection of the DIG-labeled riboprobe was done using an 
alkaline phosphatase-linked anti-DIG antibody (Boehringer Mannheim). All 
subsequent steps were carried out at room temperature. The slides were washed in 
detection buffer (DB; 100 mM Tris pH 7.5, 150 mM NaCl, 0.1% Tween), incubated 

25 with SuperBlock buffer (Biogenex) for 10 minutes and then incubated for 2 hours 
with the antibody (1:75 dilution) in DB and then washed three times in DB. The 
chromagen solution OFast Red, Sigma biochemicals) was then added and the slides 
were left to incubate for 30 minutes. The reaction was stopped by washing in 10 mM 
Tris pH 8.0, 1 mM EDTA. The cells were mounted using SlowFade (Molecular 

30 probes) and examined on a fluorescent microscope connected to a CCD camera. 

Microscopy 

Purified eosinophils were incubated in RPMI 1640 media 
supplemented with 0.5% fetal bovine serum in the presence of the compound to be 
35 tested for 15 minutes in a 24-well dish. Light microscopy was performed with an 
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inverted Axio vert 25 (Zeiss) and images were obtained with a 35 mm SLR camera 
(ARIA CONTAX, Kyocera corporation) using Kodak Elite Chrome 160T film. 

Eosinophil Chemokinesis 
5 Purified eosinophils were resuspended at 3.0 X 10^ cells per ml in 

RPMI 1640 medium supplemented with 0.5% (v/v) fetal bovine serum. Compounds 
to be tested were added from lOOOX concentrated stock solutions to 100 \i\ of cells in 
a 1.5 ml centrifuge tube and incubated at room temperature for 5 minutes. 100 |jl1 of 
treated cells were then added to the top half of a chemotactic chamber (6.5mm 

10 Trans well, 3.0 |xm polycarbonate membrane from Costar) and 600 |il of RPMI, 

supplemented with 0.5% (v/v) fetal bovine serum, was added to the bottom chamber. 
After a 2 hour incubation at 37^C in a CO2 chamber, the top chamber was discarded 
and the number of cells that had migrated to the lower chamber was evaluated by 
counting the cells using an hematocytometer. For each condition tested, the number 

15 of migrating eosinophils in two chemotactic chambers was averaged. 

Eosinophil Degranulation 

Purified eosinophils were resuspended at 3.0 x 10^ cells per ml in 
RPMI 1640 medium supplemented with 0.5% (v/v) fetal bovine serum. Compounds 

20 to be tested were diluted 1:1000 to their final concentration in 300 |il of cells in a 1.5 
ml centrifuge tube. The cells were immediately transferred to a 24-well plate. After 
an 1 hour incubation at 37 °C in a CO2 chamber, the cells were removed by 
centrifugation (4 ^C, 300xg, 10 minutes). Eosinophil cationic protein (ECP) in the 
supernatant was quantified by a double antibody radioinmiunoassay (Pharmacia) 

25 following the manufacturer's protocol. 

In Vitro Eosinophil Apoptosis Assay 

Purified eosinophils were resuspended at 2.0 x 10^ cells per ml in 
RPMI 1640 medium supplemented with 10% (v/v) fetal bovine serum, 2 mM 

30 glutamine, and 100 units of penicillin and streptomycin. Compounds to be tested 
were diluted 1 : 1000 to their final concentration and the cells were incubated at 37 °C 
in a CO2 chamber for 36 hours. The extent of apoptosis in the eosinophil population 
was evaluated using the TACS™ Annexin-V-FTTC apoptosis detection kit (R&D 
systems). Non-apoptotic cells are not stained with either Annexin-V FTTC or 

35 propidium iodide. Early apoptotic cells are stained with Annexin-V FTTC but not 
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propidium iodide (green fluorescence). Late apoptotic cells are stained with both 
Annexin-V FTTC and propidium iodide (dual green and red fluorescence). Necrotic 
cells are only stained with propidium iodide (red fluoresecence). Labeled eosinophils 
were analyzed in a FACS Calibur system from Becton Dickinson. 

5 ' ■ ■ 

Example 2 : Eosinophil Expression of DP and CRTH2 

To establish which PGD2-binding receptors are expressed by human 

eosinophils, RT-PCR was performed on total RNA from human eosinophil (>95% 
purity). The identity of the PGR products was confirmed by Southern blot detection 

10 using DP receptor and CRTH2-specific probes. 

CRTH2 mRNA was detected in eosinophils from four donors while 
DP mRNA was detected in only two of the four donors (Figure 1). The identity of the 
cell type expressing DP receptor as an eosinophil was confirmed by in situ 
hybridization. DP antisense hybridized only to cells showing the characteristic bi- 

15 lobal nucleus of eosinophils. 

Exam ple 3: PGP? Induced a Change in Eosinophil Morpholo gy Through GRTH2 
PGD2 (<10 nM) induced dramatic changes in cell morjphology within 

minutes. Vehicle-treated eosinophils were spherical and only weakly adhered to the 
20 culture dish. In contrast, eosinophils treated with PGD2 become flat, assumed an 
Amoeba-like shape and showed round structures which may represent secretory 
vesicles (panel 2b). This effect was seen on the majority of eosinophils from all 
donors analysed (n=20). PGD2-treated eosinophils reverted to a spherical shape 

within six hours. These cells were resistant to morphology changes after a second 
25 PGD2 challenge (data not shown). 

The DP receptor selective agonist, BW245C, at concentrations as high 
as 1 \xM did not affect eosinophil shape (panel 2c). In contrast, a CRTH2 selective 
agonist, DK- PGD2, induced a morphological change identical to that observed with 
PGD2 (panel 2d). Known activators of eosinophils such as platelet activating factor 

30 (PAF) (panel 2e) as well as interieukin-5 (11-5) (panel 2f) also lead to a rapid change 
in eosinophil morphology. Other prostanoid receptor agonists such as PGE2 (EP 

receptors) as well as fluprostenol and latanoprost (FP receptor) did not cause, even at 
IjM doses, aiiy alterations of eosinophil morphology (data not shown). These data 
sugjgest that the morphological changes induced by PGD2 and DK- PGD2 on 
35 eosinophils are mediated through the CRTH2 and not the DP receptor. 
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Example 4: PGD9 Increases Eosinophil Chemokinesis Through CRTH2 

PGD2 increased cell motility in the absence of a cheinical gradient, a 
process defined as chemokinesis. PGD2 was not observed to exert a chemotatic 
5 effect. Overall, the data indicates that PGD2 modulates eosinophil chemokinesis in a 

DP-independent mariner and most likely through the CRTH2. 

Chemokinesis was measured by incubating eosinophils with PGD2 for 

5 minutes prior to their loading in the upper chamber of a chemotactic unit lacking a 
chemoattractant in the lower chamber. PGD2 at a concentration of 100 mM increased 

10 eosinophil chemokinesis by 6-fold compared to cells treated with vehicle only (Figure 
3). PGD2 at concentrations of 10 nM and 1 pM increased eosinophil chemokinesis 

by a factor of 5 and 9-fold respectively (data not shown). PAF and 11-5 were also able 
to stimulate eosinophil chemokinesis (Figure 3). (See, Wardlaw, et al, /. Clin. Invest. 
78:1701, 1986, Schweizer, et a/., 7. Leukoc. Biol 59:347, 1996.) DK- PGD2 but not 

15 BW245C was effective in stimulating eosinophil chemokinesis (Figure 3). EP and FP 
receptor agonists, PGE2, fluprostenol and latanoprost failed to modulate eosinophil 

migration. 

Chemotaxis was measured by adding PGD2 to the bottom chamber of 

a chemotactic unit containing eosinophils in the top chamber. In contrast to 
20 chemoattractants such as PAF and eotaxin, PGD2 was not a chemoattractant since it 

did not attract eosinophils to the lower chamber of the chemotactic unit (data not 
shown). Eosinophils pre-incubated with PGD2 (up to 1 |iM for 5 minutes to 18. 

hours) did not haye an altered chemotactic response to either PAF or eotaxin (data not 

shown). 

25 

Example 5: PGD9 Triggers Eosinophil Degranulation Through CRTH2 

Degranulation of eosinophils was assayed by challenging freshly 
isolated eosinophils with PGD2 and measuring release of the eosinophil cationic 
protein (ECP) into the media using an ECP-specific radioinmiunoassay. PGD2 at a 
30 concentration of 10 to 100 nM significantly increased the release of ECP from 
eosinophils into the media (Figure 4). 

Among the donors, ^ broad range in PAF- and PGD2-induced release 
of ECP is seen. On average, the extent of PGD2-induced ECP release is about half 
that seen with PAF. In general the extent of degranulation induced by PGD2 and PAF 
35 paralleled each other. ECP was not released as a result of necrosis as lactate 
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dehydrogenase (LDH), a marker for necrotic cell lysis, was not detected in the media 
(data not shown). We also observed the release from eosinophils of another granular 
protein, EDN, after PGD2 challenge (data not shown). DK-PGD2, but not BW245C, 

significantly increased the release of ECP (Figurp 4). FP and EP receptor agonists did 
5 not induce ECP release. As seen with eosinophil morphology changes and 
chemokinesis, PGD2 induces eosinophil degranulation by a CRTH2-dependent 

mechanism. 

Example 6 : A Selective DP Agonist Delays The Onset Of Apoptosis 
10 The ability of PGD2 to modulate apoptosis in eosinophils was 

measured by quantifying the capacity of Annexin V to bind to phosphatidylserine on 
the outer membrane of apoptotic cells. (Koopman, et a/., Blood 5^:1415, 1994.) 
Necrotic cell death was determined by propidium iodide uptake. (Darzynkiewicz, et 
aLy Cytometry 13:195, 1992.) Annexin V-FTTC and propidium iodide staining of 
15 eosinophils was evaluated by FACS analysis. 

Isolated eosinophils become apoptotic after approximately 12 hours 
when cultured in RPMI-1640 supplemented with 10% fetal bovine serum. After 48 
hours almost all eosinophils were dead (data not shown). Addition of 11-5 or PGE2 to 

the media increased the percentage of non-apoptotic eosinophils at 36 hours in culture 
20 (Figures). 

PGD2 was a weak inhibitor of apoptotic cell death while DK- PGD2 
had no significant effect (Figure 5). In contrast, the DPrspecific agonist BW245C 
significantly increased the percentage of non-apoptotic eosinophil by 17%. The 
effects of FP agonists, fluprostenol and latanoprost were not significant. 

25 

Other embodiments are within the following claims. While several 
embodiments have been shown and described, various modifications may be made 
without departing from the spirit and scope of the present invention. 
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WHAT IS CLAIMED IS: 

1. A method of assaying the ability of a test compound to 
modulate prostaglandin D2 receptor activity comprising the steps of: 
5 a) providing said test compound to an eosinophil; and 

b) measuring the effect of said test compound on either apoptosis 
or degranulation as a measure of the ability of said test compound to modulate 
prostaglandin D2 receptor activity. 

10 2. The method of claim 1 , wherein prior to said step (a) said test 

compound has been identified ais able to bind to the DP receptor and said method 
mieasures apoptosis. 

3. The method of claim 1, wherein said method further comprises 
15 providing a DP receptor agonist to said eosinophil and measuring the ability of said 

test compound to affect apoptosis. 

4. The method of claim 1 , wherein prior to said step (a) said test 
compound has been identified as able to bind to CRTH2 and said method measures 

20 degranulation. 

5. A method of assaying the ability of a test compound to 
modulate CRTH2 receptor activity comprising the steps of: 

a) identifying a conipound that binds to a human CRTH2; 
25 b) providing said test compound to an eosinophil; and 

c) measuring either eosinophil morphology, chemokinesis under 
conditions distinguishing chemokinesis from chemotactic ability, or degranulation. 

6. The method of claim 5, wherein said step (c) measures 
30 eosinophil morphology. 

7. The method of claim 5, wherein said step (c) measures 

chemokinesis. 



17 



wo 03/008978 



PCT/CA02/01112 



- - - 8. " The method of claim Srwherein said Step (c) mea^ 
degranulation. 

9. A method of assaying the ability of a test compound to 

5 modulate CRTH2 activity comprising the steps of : 

a) providing said test compound and an CRTH2 agonist to an 

eosinophil; and 

b) measuring either eosinophil morphology, chemokinesis under 
conditions distinguishing chemokinesis from chemotactic ability, or degranulation. 

10 

10. Themethod pf claim 9, wherein said step (c) measures 
eosinophil morphology. 

11. The method of claim 9, wherein said step (c) measures 

15 chemokinesis. 

12. The method of claim 9. wherein said step (c) measures 

degranulation. 

20 13. Themethodof claim 12, wherein said agonist distinguishes 

CRTH2 from the DP receptor. 

14. The method of claim 13, wherein said agonist is 13-14-dihydro- 
15-keto-PGD2. 

25 ■ 
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SEQUENCE LISTING 



<110> 



Merck Frosst Canada & Co. 



<120> 



EOSINOPHIL PROSTAGLANDIN D2 RECEPTOR ASSAYS 



<130> 



8426-1551PCT KPM/ch 



<150> 
<151> 



60/306,357 
2001-07-18 



<160> 



4 



<170> 



FastSEQ for Windows Version 4,0 



<210> 
<211> 
<212> 
.<213> 



1 

20 
DNA 

Artificial Sequence 



<220> 
<223> 



Oligonucleotide primer 



<400> 1 

acaactcgtt gtgccaagcc 20 

<210> 2 
<211> 20 . . 
<212> DNA • 

<213> Artificial Sequence 
<220> 

<223> Oligonucleotide primer 
<400> 2 

gcatcgcata gaggttgcgc . . 20 

<210> 3 
<211> 22 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Oligonucleotide primer 



<210> 4 . 
<211> 20 . 
<212> DNA 

<213> Artificial Sequence 
<220> 

<223> Oligonucleotide primer 
<400> 4 

caggtgagca cgtagagcac * 20 



<400> 3 

ctacaatgtg ctgctcctga ac 
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7TM receptor mathematically manipulating the physicochemical descriptoi^ of step (iv) to obtain a simplified measure of the physic- 
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physicochemical descriptor or. if relevant, the simplified measure for the first 7TM receptor with the physicochemical descriptors 
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A pseudo-sequence method for comparing 7TM receptors with respect to the 
physicochemical properties of their binding sites 

Field of the invention 

5 The present invention relates to comparisons of 7TM receptor proteins with respect to 
the physicochemical properties of amino acid residues in the binding site without 
information from ligands interacting with the receptor. 

Baclcground 

10 7TM receptor proteins constitute the largest family of biological targets for current 
drugs. Drug development of ligands for these receptors has not been able to benefit 
from efficient structure-based drug design technologies due to the lack of detailed 
structural information. Determination of the exact structure of the binding site of 
transmembrane receptors is not generally amenable to X-ray crystallography or NMR 

15 as for soluble proteins, due to the inherent difficulties in producing large amounts of 
active membrane protein and in the case of X-ray crystallography producing crystals of 
these. Accordingly, the 3D structure of bovine rhodopsin is the only receptor that so far 
has been solved by X-ray crystajlography at atorhic resolution (Palczewski. K. et al, 
Sc/ence, (2000) 289. 739-745). 

20 

Despite the low sequence similarities, homology models of other 7TM receptors of the 
Class A type have been derived based on the bovine rhodopsin structure. Identification 
of binding sites and binding modes of ligands to 7TM receptors are supported by 
combining infonmation from other experimental techniques; for example, site-directed 

25 mutagenesis, metal-ion site engineering, substituted cysteine accessibility method, site 
directed spin labelling and phptoaffinity labelling. For references on designing doigs 
for 7TM receptors or G-protein coupled receptors (GPCRs), see: Drug design 
strategies for targeting G-protein coupled receptors. Klabunde.T.; Messier, G. 
ChemBioChem (2002) 3, 928-944; G-Protein coupled receptors: Models, mutagenesis, 

30 and drug design. Bikker, J.A.; Trumpp-Kallmeyer, S.; Humblet, C. Journal of Medicinal 
Chemistry (1 998) 41 , 291 1 -2927; Modelling G-protein coupled receptors for dmg 
design. Flower, D.R. Biocimica et Biophysica Acta (1999) 1422, 207-234; Locating 
ligand-binding sites in 7TM receptors by protein engineering, T.W. Schwartz. Cum 
Opin. Biotechnol. (1 994) 5. 434-444. 
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Calculated similarity scores to infer the overall "global" protein sequence homology are . 
known (Needleman S. & Wunsch C. L. J. MoL Biol. 48, 444-453 1970; Smith T.F. & 
Watennan M. J. MoL Biol. (1981) 147, 1 95-95;. Pearson W;R. & Lipman D.J. Proa Natl, 
Acad. ScL USA, (1 988) 85, 2444-48 ; Altschul S.F. J. MoL BioL (1 991 ) 219, 555-565; 
5 Altschul S.F. J. MoL EvoL (1993) 36, 290-300). Traditional evolutionary phylogenetic 
analyses considering the protein* "global" sequences are also known. Conventional 
sequence similarity scores and phylogenetic analysis of protein sequiences are based 
upon statistics on how frequently different amino acids are evolutionary replaced with 
other amino acids as reported in substitution matrices such as the BLOSUM (derived 
10 form blocks of aligned sequences) and PARM (point accepted mutations) series 
(Henikoff S. & Henicoff J.G. Proc. NatL Acad. Sci USA, (1992) 89, 10915-919; 
Henikoff S. & Henicoff J,G. Profe/ns (1993) 17, 49-61 ; Dayhoff M. Schwartz R.M. & 
Orcutt B.C. Atlas of protein Sequence and Structure (1978) 5, 345-52) and they reflect 
how nature has happened to replace them during the evolutionary period. 

15 

As illustrated in the conventional phylogenetic analysis of the GPR44 (CRTH2) 
receptor (Figure 1 ), the following relationships with a number of receptors are revealed 
according to reference J. Exp, Med. Volume (2001) 193, 255-261 .Notably, the ATI 
and AT2 receptors are not Identified according to this evolutionary relationship model. 
20 We will later show how different the relationships of r^ecieptors to GPR44 with respect to 
the ligand-binding properties in the binding site are. 

Such evolutionary-based similarity scores of sequences are different from similarity 
scores based on actual physicochemical properties associated with the individual 
25 amino acids. Orily in the cases where phylogenetically very closely related proteins are 
compared, will the analyses yield a similar result. 

Other means to compare mono-amine related 7TM receptors based on 
chemogenomics input have been devised (A novel chemogenomics knowledge-based 
30 ligand design strategy - Application to G-protein-coupled receptors, E. Jacoby. Quant. 
Struct'Act Relet. {20Q^) 20, ^^5-^23). 

Sequence similarity (as defined in various substitution matrices) versus chemical 
sequence similarity will usually not produce comparable similarity scores or identify the 
35 same related receptors. In this invention, selected amino acid residues defined as 
being part of the 7TM receptor binding site and constituting an amino acid pseudo- 
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sequence are assigned physiochemicar descriptors, which are compared and ranked 
according to a given receptor of interest 

There exists a need for a novel method for comparing 7TM receptors with respect to 
5 the physicpchemical properties of thieir binding sites, without using information from a 
ligand. This is especially important when considering 7TM receptors for which nd 
ligands are l<nown e.g. orphan receptors. 

Summary of the invention 

10 An understanding of the physicochemical properties of a binding site will assist in 
designing ligands that may bind to a receptor. The present invention describes 
methods for comparing and/or ranking 7TM receptors according to the physicochemical 
properties of their binding sites, allowing similarities and differences between said 7TM 
receptors to be identified. 

15' 

In one embodiment, the method according to the present invention relates to a pseudo- 
sequence method for comparing a first TTiVI receptor with one or more further 7TM 
receptors with respect to the physicochemical properties of their binding sites, the 
method comprisinjg the steps of: 

20 

i) optionally, aligning part of or all of the amino acid sequence of the first 7TM receptor 
vvith part of or all of the amino acid sequence of the one or more further TT M receptors, 

ii) selecting, in a sequential or non-sequential order, at the most 12 amino acid 

25 residues per helix, and at ttie most 12 amino acid residues in one or more extracellular 
loops, which are involved in one or more binding sites of each 7TM receptor, 

ill) forming a pseudo-sequence comprising at the most 50 amino acids from the 
selected sequential or non-sequential amino acid residues, 

30 

iv) for each 7TM receptor assigning one or more physicochemical descriptors to the 
amino acid residues of the selected amino acid pseudo-sequence involved in one or 
more binding sites, 
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4 

v) optionally, for each 7TM receptor mathematically manipulating the physicochemical 
descriptors of step iv) to obtain a simplified measure of the physicochemical properties 
of the binding site, 

5 vl) for each 7TM receptor generating a similarity score as defined herein by conriparing 
the physicochemical descriptor or. If relevant, the simplified measure for the first 7TM 
receptor with the physicochemical descriptors or. If relevant, the simplified measures 
for the one or further 7TM receptors, 

10 vii) optionally, ranking the 7TM receptors with respect to the physlcocherhlcal 
properties of their binding sites according to the similarity scores obtained In step vi). 

In a second embodiment, the present invention describes a method for classifying 7TM 
receptors according to the physicochemical properties of their binding sites 

15 

The invention further relates to dmg discovery methods for identifying Uganda, which 
bjnd to a first 7TM receptor and bind or potentially bind to one or more further 7TM 
receptors. Additionally, methods are described for Identifying a lead compound or a 
potential lead compound for a 7TM receptor. Furthermore, the present Invention 
20 describes a drug discovery method for constructing a phannacophore model for a 7TM 
. receptor. 

Importantly, the methods of the present Invention can be carried out using only 
infonfnatloh from the 7TM receptor, I.e. without any infomriatlon on a ligand or oh llgand- 
25 receptor interactions. 

Detailed description of the invention 

The following abbreviations and definitions are used throughout the description. 

30 A glossary of ternis used in Medicinal Chemistry or Computational Drug Design is 
found in Annual Reports in Medicinal Chemistry, Volume 33, Division of Medicinal 
. Chemistry of ACS, Academic Press, pages 385-409. For a discussion of the difference 
between "Leads" and "Drugs", see T.I. Oprea et al, J. Chem. Inf. Comput. Sci., 2001 , 
41,1308-1315. 

35 
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7TM receptor. I.e. a 7-transmembranal (7TM) receptor having seven a-heiices that 
span the cell membrane and are usually coupled to G-protelns. A list of specific 7TM 
receptors is given herein. The term is used interchangeably with "receptor . and one 
skilled in the art will understand the meaning according to the context. 

A "binding site" is a region of a biological molecule (e.g. a protein such as a TTM 
receptor) to which a iigand may bind. A binding site comprises one or more amino acid 
residues arranged in a particular geometry so as to proyide an environment with a 
specific arrangement of charged, polar or non-polar regions, which can interact with a 
Iigand. The binding site could represent the region encompassing the entire iigand, or 
consist of the main ligand-binding site, or be a subsite engaged in interactions with a 
part of the iigand. 

An 'amino acid pseudo-sequence" is defined as selected amino acid residues in 
sequential or nort-sequential order, involved in one of more Iigand binding sites in TTM 
receptors. In comparing pseudo-sequences for different 7TM receptors the positions in 
the pseudo-sequence correspond to the same generic numbers. E.g. a generic 
pseudo-sequence could be exempilfied with amino acids positioned in 1-04, 1-15, 1-^20, 
11-02, 11-07, 11-21, 111-06, 111-13 VII-02, VII-09. 

"Binding site modeF, I.e. a model achieyed by computer-assisted molecular design 
(CAMD) that describes how a Iigand binds to a TTM receptor. 

A "bitmap" is defined as a string of physlcochemical descriptor values used to describe 
certain chemical features of amino acid residues of interest. A bitmap may be derived 
from all physicochemical descriptors of the amino acid residues of the binding site or 
from a simplified measure of the physicochemical descriptors. 

"Ctiemogenqmics," i.e. correlation of chemical features of known biologically active 
compounds (i.e. ligands) with various biological targets. Cf. "Methodologies used for 
analysis of common drug shapes" G. W. Bemis and M. A. Murcko, J. Med. Cliem., 
1996. 39. 2887-2893. 

"Computer-assisted drug design (CADD)", i.e. computational techniques used to 
discover, design, and optimise small organic compounds that are biologically active 
compounds with a putative use as drugs. 
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"Design", i.e. application of all techniques leading to the discovery of new chemical 
entities (e.g. ligands) with specific properties such as affinity for a given receptor. 

5 T&ra/y or c/)e/77/ca////)rao^J.e. collection of liga^^^^ 

parallel synthesis or represent a collection of historic or commercial conipdunds. Cf. B. 
A. Bunin et al, Ann. Rep, Med. Chem, 1999, 34, 267-286. 

"Ligand", I.e. small organic compound that might display affinity for a biological target 
10 molecule such as a 7TM receptor. Throughout the specification ihe term ligand is 

sometimes used equivalently with the term small organic compound. A person skilled in 
the art will understand the meaning according to the context. 

"Main ligand-binding site" refer to the binding site located between TM-III, IV, -V, -VI, 
15 and VII in 7TM receptors which corresponds to where, for example, the main part of 
retinal is found in the rhodopsin structure and to where ligands have been mapped to 
bind in a variety of 7TM receptors. Residues, which line this site, are those that are 
generally involved in ligand binding, and therefore those, which are preferred residues 
in alignments in the present invention. 

20 

"Phamiacophore modeP, i.e. a model describing the combination of steric and 
electronic features of a ligand that are necessary for interaction with a specific 7TM 
receptor which may trigger or block its biological response. 

25 "Physicochemical descriptors" can be experimentally derived and/or theoretically 
calculated. The descriptors reflect 7TM receptor-ligand interaction features of the 
amino acid residues, i.e. they may reflect hydrophobic properties, electronic properties, 
steric properties or hydrogen bonding capabilities and other properties of importance 
for ligand-protein interactions. Some descriptors can be seen to reflect combinations of 

30 such properties, especially combinations of electronic and steric features. The 
descriptors also include dummy parameters or indicator variables, e.g. 1 and 0, to 
denote absence or presence of certain properties such as absence or presence of 
aromatic side chains, hydrophobic side chains, negatively charged side chains, 
positively charged side chains, polar side chains, hydrogen-bond donating side chains, 

35 hydrogen-bond accepting side chains or other selected features 
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"Receptor modet, i.e. a 3-dimensional model of a biological target molecule such as a 
7TM receptor based on Information from structurally known analogous proteins 
(homology model) and complementary data such as structure-activity data of ligands or 
antibodies binding to the biological target molecule and mutational studies. 

5 

A "similarity score" is a mathematical indicator of the similarity between two (potential) 
binding sites. Expressions of similarity scores include the Tanirhoto coefficient, the 
Tversky coefficient and the Euclidian distance measure. TTM receptors having close 
similarity scores have a high degree of similarity of their (potential) binding sites. 

10 

A "small organic compound" is intended Jo indicate a small organic molecule of low 
molecular weight such as below 1000. The small organic compounds of specific 
interest in the present context are those that are capable of interacting with a 
membrane-associated protein such as a 7TM receptor, in such a way as to modify the 
15 biological activity thereof . 

"Structure-based drug desiign", i.e. using protein structural information from e.g. X-ray 
crystallography or NMR spectroscopy to assist in the design of therapeutic compounds 
- mostly of inhibitors which bind to enzymes. Cf. M.A. Murcko at al, Ann. Rep. Med. 
20 C/7em, 1999, 34, 297-306. 

In the method according to the present invention, the amino acid residues of a 7TM 
receptor are assigned physicochemical descriptors. It is then possible to compare 
_ (different transmembrane proteins such as 7TM receptors with respect to the 
25 physicochemical properties of their binding sites. The process of comparing biological . 
target proteins based on the physicochemical properties of their binding site will be 
referred to "physicogenomics" herein. The following steps are involved in 
physicogenomics: 

• Sequence search 
30 • Alignments 

• Analysis of binding site residues important for ligand binding and recognition 

• Comparison of physicochemical properties of binding site 

In other words, the method according to the present invention relates to a method for 
35 comparing a first 7TM receptor with one or more further 7TM receptors with respect to 
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the physicochemical properties of their binding sites, the method comprising the steps 

i) optionally, aligning part of or all of the amino acid sequence of the first 7TM receptor 
5 with part of or ail of the amino acid sequence of the one or more further 7tM receptors, 

ii) selecting, in a sequential or non-sequential order, at the most 12 amino acid 
residues per helix, and at the most 12 amino acids in one or mote extracellular loops, 
which are involved in one or more binding sites of each 7TM receptor, 

10 

iii) forming a pseudp-sequence comprising at the most 50 amino acid residues from the 
selected sequential or non-sequential amino acid residues, 

iv) for each 7TM receptor assigning one or more physicochemical descriptors to the 
15 amino acid residues of the selected amino acid pseudo-sequence involved in one or 

more binding sites, 

v) optionally, for each 7TM receptor mathematically manipulating the physicochemical 
descriptors of step iv) to obtain a simplified measure of the physicochemical properties 

20 of the binding site, 

vi) for each 7TM receptor generating a similarity score as defined herein by comparing 
the physicochemical descriptor or, if relevant, the simplified measure for the first 7TM 
receptor with the physicochemical descriptors or, if relevant, the simplified nfieasures 

25 for the one or further 7TM receptors, 

vii) optionally, ranking the 7TM receptors with respect to the physicochemical 
properties of their binding sites according to the similarity scores obtained in step vi). 

30 Importantly, this comparison or ranl<ing can be carried out using only information from 
the 7TM receptor, i.e. without any information on a ligand. Therefore, the present 
invention also relates to a method as described herein, wherein the comparison is 
made without using data related to binding affinity of a ligand to a 7TM receptor. 



35 



The present invention may also be used to classify 7TM receptors according to the 
physicochemical properties of their binding sites. Hence, the present invention 
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describes a method for classifying 7tM receptors according to the physicochemical 
properties of their binding sites. This classification may also be earned out using only 
information from the 7TM receptor, i.e. without any infomiation on a ligand. Hence, the 
classificatioh may be made without using data related to binding affinity of a ligand to a 
5 7TM receptor. 

Method for obtaining lead structures fornew 7TM receptors 

It is generally recognized that structurally related small organic molecules often bind 
the same biological target proteins. Cf. "Do Structurally Similar Molecules Have Similar 
10 Biological Activity?" Yvonne C. Martin, James L. Kofron, and Linda M. Traphagen, J. 
Med. Chem. (2002) 45, 4350-4358. Conversely, related biological targets often bind the 
same or similar small organic molecule llgands. Notably, there are exceptions from 
these general rules - othenwise development of receptor-selective ligands would have 
been an impossible task. 

15 

Having identified which 7TM receptors resemble each other with respect to the 
physiochemical environment in the binding site as described herein, it is possible to 
utilise known ligands which interact with one receptor as chemical starting points (so 
called chemicalleads) for drug development on related receptors. 

20 ' ■ ' 

in other words, the present Invention also relates to a drug discovery method for 
identifying ligands, which bind to a first 7TM receptor and potentially bind to one or 
more further 7TM receptors, the method comprising the steps of 1) to vii) as defined 
above and the further steps of 

25 

viii) selecting from one to about 100 further 7TM receptors which have the closest 
similarity scores to the first 7TM receptor, 

ix) identifying ligands which potentially bind to those further 7TM receptors selected in 
30 step viii) by selecting ligands that bind to the first 7TM receptor. 

In addition, the method described above may be combined with a screening step, to 
detenmine ligands, which do bind to a 7TM receptor. Therefore, the present invention 
additionally relates to a drug discovery method for identifying ligands which bind to a 
35 first 7TM receptor and to one or more further 7TM receptors, this method comprising 
the steps of 1) to vii) as defined above and the further steps of: 
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viii) selecting from one to about 100 further 7TM receptors which have the closest 
similarity scores to the first 7TM receptor, 

5 ix) screening ligands that bind to the first 7TM receptor against the selected 7TM 
receptors of step viii). 

The process of transfer of a chemical starting point from one protein target to another 
related protein target is often referred to as chempgenomics. Thus, an efficient 

10 physicogenomic method to compare 7TM receptors having known ligands (known or 
potential drug molecules) with novel receptors lacking identified ligands:allows for 
possibilities to identify lead structures for drug development since no previous 
information regarding ligands binding to the new receptor under investigation is 
needed. Furthermore, comparison of 7TM receptors having ligands with orphan 

15 receptors lacking identified endogenous ligands allows for the identification of lead 
structures for drug development also on orphan receptors since no previous 
information regarding ligands binding to the receptor under investigation is needed. 
Consequently, known ligands of closely related receptors could serve as good 
chemical starting points to identify lead structures against a receptor for which no 

20 agonist or antagonist are known. 

In a furiiher aspect, the present invention relates to a drug discovery method for 
identifying a potential lead compound for a first 7TM receptor, the method comprising 
the steps of i) to vii) as defined above and the further steps of 

25 

viii) selecting from one to about 100 further 7TM receptors which have the closest 
similarity scores to the first 7TIVI receptor, 

ix) identifying ligands that bind to said one or more further 7TM receptors to construct a 
30 library including a potential lead compound for the first 7TM receptor. 

The method for identifying a potential lead compound may additionally be linked to a 
screening step, so that libraries containing potential lead compounds may be narrowed 
down to give lead compounds. Hence, the present invention further relates to a drug 
35 discovery method for identifying a lead compound for a first 7TM receptor, the method 
comprising the steps of I) to vii) as defined in herein and the further steps of 
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viii) selecting from one to about 100 further 7TM receptors which have the closest 
similarity scores to the first 7TM receptor, 

5 ix) identifying ligands that bind to said one or more furthier 7tM receptors to construct a 
library, and 

x) screening said library against the first 7TM receptor to identify a lead compound for 
the first 7TM receptor. 

10 

Pharmacophore model ^. 

Having established physicogenomic relationships between 7TM receptors one can also 
derive pharmacophore models based on ligands, which are known, to bind to 7TM 
receptors related to the receptor under investigation. This allows for a more cost- 
15 effective process of retrieving and screening a limited number of compounds than 

under a conventional high-throughput screening (HTS) campaign usually conducted in 
larger pharmaceutical Industries when looking for lead structures for a novel biological 
target, 

20 Furthermore, a pharnniacophore model that can be used for in silico screening or design 
of focused chemical libraries could be derived from analogously identified structures. 

Hence, the present invention relates to a drug discovery method for constructing a 
. .. phannacophore model for a first 7TM receptor, the method comprising the steps of i) to 
25 vii) as defined in herein and the further steps of 

viii) selecting from one to about 1 00 further 7TM receptors which have the closest 
similarity scores to the first 7TM receptor, 

30 ix) identifying ligands that bind to said one or more further 7TM receptors to construct a 
pharmacophore model. 



35 



In one embodiment of the present invention the first 7TM receptor is one for which no 
ligands have been identified. Additionally, the first 7TM receptor may be an orphan 
receptor. 
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7TM receptors, which have similarity scores closest to each other, have most similar 
physicochemlcal properties of their binding sites. Therefore, when selecting 7TM 
receptors based on their similarity score, with the aim of choosing those with similar 
physicochemical properties of their binding sites, it is important to select those with the 
5 closest similarity scores. The number of further 7TM receptors which are selected in 
step vii) (above) may be from one to 50, from one to 25 or from one to 15. 

Having defined the relevant binding site and assigned the proper physicochemical 
descriptors to selected amino acid residues involved in the binding site, the 7TM 

10 receptors are compared to each other by a suitable computerised mathematical model, 
which is capable of comparing-a large number of receptors by an effective algorithm. 
Preferably, all identified 7TM receptors should be aligned and compared based on the 
physicochemical descriptors of pseudo-sequences derived from selected amino acid 
residues involved in the binding site. Accordingly, the present invention relates to any 

15 of the methods as described herein, wherein the method is executed by a computer 
under the control of a program and the computer includes a memory for storing said 
program. 

Generic numbering system for 7TM receptors 

20 The 7TM receptor superfamily is composed of hiahy hundreds of receptors that may be 
further divided into smaller sub-families of receptors. The largest of these sub-families 
of 7TM receptors is composed of the rhodopsin-like receptors (also termed the family A 
receptors), which are named after the light-sensing molecule from our eye. The 
receptors are integral membrane proteins characterized by seven transmembrane 

25 (7TM) segments traversing the membrane in an antiparallel way, with the N-terminal on 
the extracellular side of the membrane and the C-terminal on the intracellular side. 
Within the membrane- embedded part and in some cases in the membrane proximal 
parts embedded in the aqueous environment surrounding the cell membrane, the 
polypeptide adopts a helical secondary structure. The lengths, and the beginning, 

30 centre and ends relative to the lipid bilayer membrane of these helices may be 

deduced from solved three-dimensional structures of the receptor proteins (Palczewski 
K. et al., Science, (2000) 289 (5480), 739-45). However, since the three-dimensional 
structure of only a single receptor has been solved to date, the helical lengths, and the 
beginning, centre and ends relative to the lipid bilayer membrane of each of the seven 

35 helices may be dissected by sequence analysis (J.M. Baldwin, EMBO J. (1993) 12(4), 
1693-703; J.M. Baldwin et al., J. Mo/ 8/o/, (1997) 272(1), 144-64). 
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A useful tool in the analysis of 7TM receptors is the generic numbering system for. 
residues of 7TM receptors. Within the many hundred members of the rhodopsin-like 
receptor family, a number of residues especially within each of the transmembrane 
5 segments are highly, but not totally, conserved - termed key residues. These residues 
may be used to direct an alignment of the primary protein sequences within the 
transmembrane segments together with other standard principlies and techniques, well 
known to persons skilled in the art (for example hydrophobicity plots). Additionally a 
number of other residues occur within the transmembrane segments that are highly 
10 conserved, and these may be used to further direct an alignment of the transmembrane 
segments. These are particularly useful when a given key residue in a transmembrane 
segment has been substituted through evolution by another amino acid of a dissimilar 
physiochemical nature. 

15 However, due to differences in the length of especially the N-terminal segment, the 
third intracellular loop 3 and the C-terminal segment, residues located at presumably 
structurally corresponding positions in different 7TM receptors are numbered differently 
in different receptors. However, based on the conserved key residues in each 
transmembrane segment, a generic numbering system has been suggested (JM 

20 Baldwin, EMBO J. (1993) 12(4), 1693-1703; TW Schwartz. Curr. Opin. Biotech. (1994) 
5, 434-444). On the basis of the key residues present in the receptor family, the 
transmembrane segments are generically numbered. For example, in TM-M the highly 
conserved acidic function, aspartate (Asp) in the rhodopsin-like family is given the 
generic number 10, i.e. AsplhlO, on_the basis of its.position in the helix. All other 

25 residues in the helix are hence numbered on this basis. In Figure 2.1 a schematic 
depiction of the secondary structure of a rhodopsin-like 7TM receptor is shown with 
one or two conserved, key residues highlighted in each transmembrane segment: 
Asnl:18; Aspll:10; Cyslll:01 and Arglll;26; TrplV:10; ProV:16; ProVI:15; ProVII:17. 

30 In relation to the present invention, it is important that residues involved In, for example, 
ligand binding can be described in this generic numbering system. For example, in the 
i92-adren(Brgic receptor the main anchor point for the catecholamine agonists is a highly 
conserved aspartate 113 residue in TM-III (Asplll:08). Interestingly, the aspartate is 
located one helical turn deeper in the receptor structure to the counter-Ion (Glulll:04) of 

35 the Schiff base in rhodopsin. Two serine residues at positions 204 and 207 in TM-V 
(SerV:09 and SerV:12) were identified as contact points for the hydroxy-functionality of 
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the catechol ring and were suggested also to be part of the control of the transition 
between the inactive and active state of the receptor. Moreover phenylalanine 290 and 
asparagine 293 on TM-VI (PheVI:17 and AsnVI:20 respectively) were Identified as 
Interacting with the catechol ring and the ;9-hydroxy of the agonists respectively. In 

5 catecholamine and 5-hydroxytryptamlne receptors, residue Vll:06 has been 
consistently identified as an interaction point for partial agonists and antagonists. 
Interestingly, a y?2-selectlve non-catechol agonist TA-2005 was suggested to Interact 
with tyrosine 308 on TM-VII (TyrVII:02), located one helical turn above Vll:06, and was 
unaffected with respect to receptor activation by mutation of either serine residues 204 

10 or 207 suggesting that different agonist activation modes may exist even within a 
particular receptor. Hence the agonist-binding site in the y92-adrenergic receptor is 
located between the transmembrane segments at a general position rather similar to 
the retinal binding site in opsin. Comparing the binding site of the ;ff2-adrenergic 
receptor described above using the generic nomenclature to, for example, the 

15 dopamine D2-receptor, quickly identifies residues Asp114 in TM-III (Asplll:08), Ser194 
and Ser197 in TM-V (SerV:09 and SerV:12). Phe390 and Hls393 {PheVI:17 and 
HlsVI:20) hence suggesting that certain key residues In the binding of the 
catecholamines to the iS2-adrenergic receptor are conserved to for example the 
dopamine receptoir. SImilariy the described binding site - or other residues suspected 

20 to be irriportant in defining for example a binding site of a llgand - can easily be 
compared to even several and very distantly related receptors at the same tirrie, where 
little or no conservation of the particular binding site of the endogenous catecholamine 
bindlnig site exist, but where the corresponding residues nevertheless may be Important 
for binding, for example, artificial npn-peptlde antagonists. 

25 

It is only in the riiodopsin-like receptor family that a generic numbering system has 
been established; however, it should be noted that although the sequence identity 
between the different families of 7TM receptors is very low, it is believed that they may 
share a more-or-less common seven helical bundle structure. Hence an analogous 

30 system may be developed for the other families of 7TM receptors, for example the 
family B class of receptors, composed of, among others, the glucagon receptor, the 
glucagon-IIke peptide (GLP) receptor-1 , the corticotropin releasing factor (CRF) 
receptor-1, vasoactive intestinal peptide (VIP) receptor, pituitary adenylate cyclase- 
activatlng polypeptide (PACAP) receptor etc. (J.W. Tarns et al., Receptors Channels 

35 (1998) 5(2), 79-90). Again, on the basis of the key residues present In the family B 
class of receptors, the transmembrane segments are generically numbered (Figure 
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2.2). For example, in TM-I, the highly conserved hydroxy function, serine (Ser) is given 
the generic number 8, i.e. Serl:08 on the basis of its approximate position in the helix; 
in TM-ll the highly conserved histidine (His) Is given the generic number 6 i.e. Hisll:06; 
in TM-III the highly conserved cysteine (Cys) is given the generic number 1 i.e. 
5 Cyslll:01 ; in TM-IV the highly conserved proline (Pro) is given the generic number 13 
i.e. ProlV:1 3; in TM-V the highly conserved asparagine (Asn) is given the generic 
number 14 i.e. AsnV:i4; In TM-VI the highly conserved leucine (Leu) is given the 
generic number 9 i.e. LeuVI:09; in TM-VII the highly conserved glycine (Gly) is given 
the generic number 13 i.e. GlyVII:13. All other residues in the helices are hence 

10 numbered on this basis. Thus, all the techniques described in the present invention can 
be applied to the other families of 7TM receptors with minor modifications. This generic 
numbering system together with general knowledge of the 3D structure of the 7TM 
receptors and knowledge from systematic ligand binding site analysis makes it possible 
to predict or identify contact points for ligands based on the DNA sequence coding for 

15 the 7TM receptor (see Figure 2.1 and 2.2). 

7TM receptors of interest in the present Invention comprising 7 transmembrane 
domains include but are not restricted to G-protein coupled receptors, such as 
receptors for: acetylcholine, adenosine, norepinephrin and epinephrine, amylln, 

20 adrenbmedullin, anaphylatoxin chemotactic factor, angiotensin, apelin, bombesin 
(neuromedin), bradykinin, calcitonin, calcitonin gene related peptide, conopressin, 
corticotropin releasing factor, , calcium, cannablnoid, CC-chemokines, CXC- 
chemokines, CX3C-chemokinees, cholecystokinin, cortlcotropin-releasing factor, 
. dopamine, eicosanoids, endothelin, fMLP, GABAe. galanin, gastrin, gastric inhibitory 

25 peptide, glucagon, glucagon-like peptide I and II, glutamate, glycoprotein hormone (e.g. 
FSH, LSH, TSH, LH), gonadotropin releasing honnone, growth hormone releasing 
hormone, growth hormone releasing peptide, Ghrelin, histamine, 5-hydroxytryptamine. 
leukotriene, lysophospholipid, melanocortins, melanin concentrating homnone, 
melatonin, motilin, neuropeptide Y, neurotensin, nocioceptin, odor components , 

30 opiods, retinal, orexins, oxytocin, parathyroid hormone/parathyroid hormone-related 
peptides, pheromones, platelet-activating factor, prostanoids, secretin, somatostatins, 
tachykinins, thrombin and other proteases acting through 7TM receptor, thyrotropin- 
releasing hormone, pituitary adenylate activating peptide, vasopressin, vasoactive 
intestinal peptide and virally encoded receptors; and 7TM receptors coded for in the 

35 human genome and for which an endogenous receptor-ligand has or has not yet been 
assigned; in particular adenosin, galanin, CC-chemokines, CXC-chemokines, 
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melanocortin, bornbesin, cannabinoid, lysophospholipid, fMLP, neuropieptide Y, 
tachykinin, doparnine, histamine, 5-hydroxytryptamine, histamine, mas-proto- 
oncogene, melanin concentrating hormone receptor, Glucoserdependent Insulinotropic 
polypeptide, Glucagonriike peptide-1 receptor, glucagon receptor, acetylcholine. 

5 oxytocin, human herpes virus encoded receptors. Epstein Barr virus induced receptors, 
cytomegalovirus encoded receptors and bradykinin receptors; and 7TM proteins coded 
for in the human genome but for which no endogenous receptor-ligand has yet been 
assigned; preferably the galanin receptor type 1 , leukotriene B4 receptor, CCR1 , 
CCR2, CCR3, CCR4, CCR5. CCR6, CCR7, CCR8, CCR9. CCR10, CXCR1, CXCR2, 

10 CXCR3, CXCR4, CXCR5, CXCR6, CX3CR1 , melanocortin-1 receptor, melanocortin-^ 
receptor, melanocortin-4 receptor, melanocort|n-5 receptor, MCH-1 , MCH-2, GIF 
receptor, GLP-1 recieptor, bombesin receptor subtype 3, cannabinoid receptor 1, 
cannabinoid receptor 2, EDG-2, EDG-4, FMLP-related receptor I, FMLP-related 
receptor II, NPY Y6 receptor, NPY Y5 receptor, NPY Y4 receptor, NK-1 receptor, NK-3 

15 receptor, D1 receptor. D2 receptor (short), D2 receptor (long), D3 receptor, D4 
receptor, D5 receptor, 06 receptor, D7 receptor, D7 receptor, D8 receptor, duffy 
antigen; US27, US28, UL33 and U78 from human cytomegalovirus; U12 and U51 from 
human herpes virus 6 or 7, ORF74 from human herpes virus 8, and histamine HI 
receptor, MAS proto-oncogene, muscarinic Ml receptor, muscarinic M2 receptor, 

20 muscarinic MS receptor, muscarinic M5 receptor, oxytocin receptor, XCR1 receptor, 
EBI2 receptor, RDC1 receptor. GPR12 receptor or GPR3 receptor, and 7TM receptors 
coded for in the human genome and for which an endogenous receptor-ligand has or 
has not yet been assigned. These 7TM receptors may be studied in a monomeric or in 
— a dimeric fomi. which may be either homo-dimeric or heterodlmeric. 

25 

Aligning 7TM receptor proteins 

As appears from the above, it is important to align at least a part of the amino acid 
sequences of the receptors to be compared. In some cases this information may 
already be available and, therefore, such an alignment step may be omitted, but in 
30 other cases the method of the invention includes a step of alignment. 

Firstly, sequence databases, such as SWISSPROT, SPTREMBL, EMBL, PIR etc. are 
searched for human GPCR sequences using the Sequence Retrieval System (network 
browser for databanks in molecular biology) SRS. The Identified sequences are then 
35 aligned using conventional alignments algorithms such as ClustalW (Thompson J.D. 
Higgins D.G. & Gibson T. J. Nucleic Acids Research. (1 994) 22, 4673-80) The 
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resulting alignment is nnanually inspected and refined if necessary, so that conserved 
generic sequence signatures within the seven transmembrane 7TM helices are 
satisfied (Palczewski. K. et al, Science. (2000) 289, 739-745). 

5 Secondly, the helices of the 7TM receptor are identified based on hydrophobicity plots, 
the conserved residues within the siub-family of receptors and - for family A receptors - 
the sequence alignment to the recent published crystal structure of rhodopsin (OPSD). 

Furthermore, the invention relates to methods as described herein, wherein step i) is 
10 included and the alignment is based on a model developed for 7TM receptors. The 
7TM receptors of the present invention may be Class A, Class B, Class C or taste 
receptors - . . . — . 

In one embodiment, the invention describes a method or methods described above, 
15 wherein step i) is included and the alignment is made with respect to transmembrane 
positioning of a-hi3lices of 7TM receptors. 

A method according to the invention involves the use of a pseudo-sequence. A 
pseudo-sequence is obtained from at the most 12 amino acid residues per 7TM helix, 

20 and al the most 12 amino acids in one or more extracellular loops, sequential or non- 
sequential, involved in one of more binding site. Furthermore, such a pseudo-sequence 
may comprise at the most 50 amino acid residues. In a specific embodiment at the 
nriost 8 such as, e.g., at the most 6 amino acid residues per 7TM helix or extracellular 
loop fomri the pseudo-sequence containing at the most 40 amino acid residues such 

25 as, e.g., at the most 30 amino acid redidues. In an embodiment of the invention only 
amino acid residues from at the most 6 such as, e.g., 5 helices are included in the 
pseudo-sequence. 

A surprinsingly good reflection of the physicochemical environment in the binding site is 
30 conveyed by the relatively few amino acids invoked in the analysis of only a relatively 
short pseudo-sequence. This is supported by the fact that ligands binding to receptors 
that have been found to be related to the given receptor also could display affinity for 
the given receptor (cf. investigation of CRTH2 and identification of the related receptor 
angiotensin II and hence a angiotensin II llgand candesartan also binding to CRTH2). 

35 
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A selection of relevant amino acid residues for the pseudo-sequence can be made by a 
person skilled in the art based on current knowledge on 7TM receptors as discussed in 
the paragraph "Generic numbering system for 7TM receptors". Examples of relevant 
literature are T. Klafaunde, G. Messier, Drug design strategies for targeting G-protein- 
5 coupled receptors. Chembiochem. 3 (2002) 928-44; D. R. Flower, Modelling G-protein- 
coupled receptors for drug design, Biochimica et Biophysica Acta 1422 (1 999) 207- 
234; J. A. Bjkker, S. Trumpp-Kallmeyer, and C. Humblet G-Protein Coupled Receptors: 
Models, Mutagenesis, and Drug Design Journal of Medicinal Chemistry 41 (1998) 
2911-2927. 

10 

Although amino acid residues known to be important for small molecule ligand 
interactions of course are of interest in construction of the pseudo-sequence, the 
present invention is not limited thereto^ 

15 ; " , ■ 

Ah example of such an alignment for a random set of selected 7TM receptor human 
sequences is shown in the table below. For this illustration, the sequences and 
identification codes for the 7TM receptors are retrieved from www.qpcrorq . The amino 
acid residues, in sequential or non-sequential order, are selected from different helices 

20 located in the binding site. For GPR44, the following amino acid residues, up to six per 
helix, from TM-III, TM-IV, TM-V, TM-VI and TM-VIII. have been selected. In an specific 
example given below, the selected residues for TM-VII correspond to VI 1-02 (Leu), VII- 
06 (Thr) and VII-09 (Ala) in the generic numbering. The rest are assigned analogously 
as described above. 

25 " ' ' ' - ^ 

TM Helix: IN IV V VI VII 

GP44_HUMAN HSFFMF NTY AKFA WYHSEA LTA 

The amino acid residues, e.g. up to six per helix, in sequential or non-sequential order, 
are selected from 111-04 to VII-09 to fomi the following 22 amino acid pseudo- 
30 sequences, which are used in the alignment and subsequent comparison. 

GP44_HUMAN HSFFMFNTYAKFAWYHSEALTA 
02T1_HUMAN QHYLVGDGLSINFLFSLYAKVT 
07C2_HUMAN QIFIGCGSETEIFVLCLYSLVT 
35 B3AR_HUMAN WTDWTVSPVSSSWFFNRAFNG 
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PE24_HUMANSTLLSLTrTAASSSLWNQDIA 

APJ_HUMAN SSIFMYLAVGSTGWYHKYMFTS 

01 E1_HUMAN QMFLGDHAHACFDVFLLYATMT 

FML2_HUMANVHIDLFLTNLHFGWYEGMAISA 
5 ACTR_HUMAN IDFVLLTGMWITWFWMTFGI 

5H4_HUMAN RTDVTTISPACSAWFFNDPWLG 

NFF2_HUMAN SGQGVAIMSTWRWLWMSDYHA 

02B2_HUMAN QLFLGSNSQHVDFLVTLYAKLG 

AG2R_HUMANASVSLYASAGKNGWHQTDVMIA 
10 05U1_HUMANQVFIASSGHKIHFRSARVFLVT 

1019_HUM NLLSRTLNLHLYEFSIGSMFLT 

C3X1_HUMANTTFFFFVAQNTNGWYNIETLEA 

5H6_HUMAN WTDVCSASPVASTWFFNQAFTG 

BRB2_HUMANVNISLYLSMNLNGWFQTDTTSA 
15 02F2_HUMAN QLSLGGNSQPTNIMFCLYIKVA 

NMBR_HUMAN IPQLVGLAESIFYWNHYRSTRS 

0XI2_HUMANQMIHSMARISLSYYMISHRVNL 

NTR2_HUMAN YYHEAYLAMINVSWYHRYCYNF 

AG2S_HUMAN ASVSLYASAGKNGWHQTDVMIA 
20 5H2A_HUMANWIDVSTISIVGSSWFFNAVLVG 

GP72_HUMANSRQYLHFSHDTFLWLNVLSYHA 

ACMI^HUMANWLDYSNLWATTAAWYNVSTWYC 

CKRA_HUMANISYSFHLAAAQVGQYSLDTLSA 

TA2R_HUMANMGMIGLLGPSLSLWLLITVLLA 
25 0YD1_HUMANQMVHYARRYGVAAYAFFHRINV 

LGR5_HUMANIGSISEKYSLLNCNVASSLKLV 

05V1_HUMAN QLFWGNSHNINFWFLVYIRVS 

FSHR_HUMAN AGTVSEAAPVCLDMISAASKVH 

DADR_HUMAN WVDISTISPASSSWFFNLPFVG 
30 GRPR_HUMAN IPQLVGLAESSFFWNHYRSSRA 

AG22_HUMANFGLTMFSSTAKNGWFHTDALIG 

02D2_HUMAN RLFLGCVSDSIAFLVFLYGKVA 

04F3_HUMAN QIIHGGHSQPLDYFPMYPHKIA 

02H3_HUMAN QIFLGTAWGQSTLVLSLYGKLA 
35 HH1R_HUMANWLDYSTWVIKTANWYFFIAHIG 

CML1 HUMAN SNLIMFLSSTRFGWYHNELLTA 
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EBI2.HUMAN TAFYTYQTLLACGFYHIHMLVM 
02GCHUMAN QMALGGSAYAIRLHSNRALNVT 
ETBR_HUMAN VPQKVGLAEDLFYWLHRKlDIA 
GPRCHUMAN VFTVLDIGATKFGWYHSELITA 
6 FML1_HUMAN IHVDLFLTVGRFGWFQAGTVSA 
NK4R_HUMAN QNPIVFLAQHVIVWYHFTAYFA 
GPR81_HUMA GLLARAGTLHMFEYSVRFLLLT 
GPRB^HUMAN TVLVFALGPLRASWFACGSTLA 

10 Table. Pseudo-sequence generic numbering alignment of selected human 7TM 
receptors. 

Selection of binding site amino acid residues 

Tine alignment of a pseudo-sequence mentioned above is used to identify binding sites 
15 or potential binding sites of tfie 7JM receptor. Such identification is necessary in order 
to enable designation of physicochemical descriptors to the amino acid residues 
involved in the (potential) binding site. 

As outlined above, certain amino acids have been identified to be frequently involved in 
20 ligand interactions. Thus, amino acid residues facing the core of the 7TM bundle and 
the general 7TM ligand binding site defined by the transmembrane helices in addition 
to residue positions in the extracellular loops determined experimentally to be important 
for ligand binding are selected. Therefore, in an embodiment of the invention, the 
binding site includes amino acid residues located in one or more extracellular loops of 
25 the 7TM receptors. In a further embodiment of the invention, the binding site includes 
amino acid residues located ih one or niore subsites of the binding site and in one or 
more extracellular loops of the 7TM receptors. 

Relevant examples on identification of residue positions of importance for ligand 
30 binding and recognition can be found in: Y Yamano, K. Ohyama, S. Chaki, D.F. Guo, T. 

Inagami, Biochem. Biophys. Res. Commun. 187 (1992) 1426-1431; S.A. Hjorth, H.T. 

Schambye, W.J. Greenlee, T.W. Schwartz, J. Biol. Chem. 269 (1994) 30953-30959; 

Y.H. Feng. K. Noda. Y. Saad, X.P. Liu, A. Husain, S.S. Karnik, J. Biol. Chem. 270 

(1995) 12846-12850; K. Noda, Y. Saad, S.S. Kamlk, J. Biol. Chem. 270 (1995) 28511- 
35 28514; H. Ji, M. Leung, Y. Zhang, K.J. Catt, K. Sandberg, J. Biol. Chem. 269 (1994) 

16533-16536; Y. Yamano, K. Ohyama, M. Kikyo, T. Sano, Y. Nakagomi, Y. Inoue, N. 
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Nakamura, I. Morishima, D.F. Guo, T. Hamakubo et al. J. Biol. Chem. 270 (1995) 
14024-14030; K; Noda. Y. Saad, A. Kinoshita, TP, Boyle. R.M. Graham, A. Husain, 
S.S. Karnik. J. B/d/. C/?em. 270 (1995) 2284-2289; HT. Scham^ 
Weinstock, T.W. Schwartz, Mol. Pharmacol. 47 (1995) 425-431); V. Nimla. W. Zheng, 
5 R. Sothinathan, K. Sandberg, Br. J. Pharmacol. 119 (1996) 15p5-1507. 

Physicochemical descriptors 

Various physiochemical descriptors have previously been employed to classify and 
describe chemical features of peptides. In a method according to the present invention 
10 physicochemical descriptors are applied to amino acid residues located in or in the 
vicinity of the (potential) binding site. Use of such descriptors enables calculation of a 
similarity score between 7TM receptors so that a comparison of the individual 7TM 
receptors can easily be made. 

15 Thus, structure-activity studies on peptides have shown the relevance of various 

physicochemical descriptors assigned to individual amino acid residues in describing 
biological properties by Quantitative Structure-Activity Relationships (QSAR), Principle 
Component Regression (PGR) and Partial Least-Squares (PLS) analysis. Such studies 
using various descriptors can be found in (among others): Simple parameterization of 

20 non-proteinogenic amino acids for. QSAR of antibacterial peptides. Lejon, Tore; 
Svendsen, John S.; Haug, Bengt E. Journal of Peptide Science (2002) 8, 302-306; 
Theory and applications of the integrated molecular transfonrn and the normalized 
molecular moment structure descriptors: QSAR and QSPR paradigms. Molnar, 
Stephen P.; King, James W. IntemaiionalJoumal of Quantum Chemistry {200^) 85, 

25 662-675; Modelling of the Amino Acid Side Chain Effects on Peptide Confomiation. 
Sak, Katrin; Karelson, Mati; Jarv, Jaak. Bioorganic Chemistry (1999) 27, 434-442; The 
application of the intermolecular force model to peptide and protein QSAR. Charton, 
. Marvin. Advances in Quantitative Structure-Property Relationships (1 999) 2, 1 77-252; 
MS-WHIM Scores for Amino Acids: A New 3D-Description for Peptide QSAR and 

30 QSPR Studies. Zaliani, A.; Gancia, E., Journal of Chemical Information and Computer 
Sciences (1999) 39, 525-533.; A validation study of molecular descriptors for the 
rational design of peptide libraries. Matter, H. Journal of Peptide Research (1998) 52, 
305-314; 3D-QSAR of angiotensin-converting enzyme inhibitors: functional group 
interaction energy descriptors for quantitative structure-activity relationships study of 

35 ACE inhibitors. Kim, Sanguk; Chi, Myung Whan; Yoon, Chang No; Sung, Ha-Chin. 
Journal of Biochemistry and Molecular Biology (1 998) 31 . 459-467; Descriptors for 
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amino acids using MolSurf parametrization. Norinder, Ulf; Svenssbn, Peter. Journal of 
Computational Chemistry (1 998) 19, 51-59; Amino Acid Side Chain Descriptors for 
Quantitative Structure-Activity Relationship Studies of Peptide Analogs. Collantes, 
Elizabeth R.; Dunn, William J., III. Journal of Medicinal Chemistry (1995) 38. 2705-13; 
5 Theoretical amino acid descriptors. Application to bradykinin potentiating peptides. 
Norinder, Ulf. Peptides (New YorKUS) (1991) 12, 1223-7; Dedicated principal 
properties for peptide QSARS: theory and applications. Skagerberg, Bert; Sjdstrbm, 
Michael; Wold, Svante. Journal of Chemometrics (1990) 4, 241-53; Multivariate 
Parametrization of 55 Coded and Non-Coded Amino Acids. Jonsson, J.; Eriksson, L; 

10 Hellberg, S.; Sjostrom, M.; Wold, S. Quant Struct-Act Relat (1989) 8. 204-209; 
Amino Acid Side Chain Parameters for Correlation Studies in biology and 
pharmacology. Fauchere, J.-L; Charton, M ; Kier, L. B.; Veriooop, A.; Pliska, V. Int. J. 
Pept Protein Res. (1988) 32, 269-278; Peptide Quantitative Structure-Activity 
Relationships, a Multivariate Approach. Hellberg, Sven; Sjostrom, Michael; 

15 Skagerberg, Bert; Wold, S\/ar\ie. Journal of Medicinal Chemistry (1987) 30, 1 126- 
1 135; Statistical Analysis of the Physical Properties of the 20 Naturally Occurring 
Amino Acids, Kidera, A.; Konishi, Y.; Oka, M.; Ooi, T.; Scheraga, J. A. J. Protein 
Chem. (1985) 4, 23-55; Relations between Chemical Structure and Biological Activity in 
Peptides.Sneath, P. H.A. J. rif}ebr.e/o/.(1966)12, 157-195. : 

20 

Such physicochemical descriptors applied to arhino acids in peptides reflect the forces 
involved in ligand-receptor interactions and, accordingly, will reflect the interacting 
properties of the amino acid side chains in proteins, especially transmembrane 
receptors such as 7TM receptors. 

25 

For example, the amino acids may be described by surface volumes and log P of side- 
chains (Norinder, Ulf; Svensson, Peter. Journal of Computational Chemistry {^998) 19, 
51 -59), ^-angles and conformational strain energies AHstrain (Sak, Katrin; Karelson, 
Mat!; Jarv, Jaak. Bioorganic Chemistry {"[999) 27, 434-442) or principle properties z 
30 (Hellberg, Sven; Sjostrom, Michael; Skagerberg, Bert; Wold, Svante. Journal of 
Medicinal Chemistry (1 987) 30, 1 1 26-1 1 35) as shown in the table below. 
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Sutface 


logP 


angle chi 


dH Strain 


Z1 Z2| Z3 






side chain 


side chain 








Ala 


A 


37.988 


-0.093 


-77.85 


0.00 


0.07 


-1.73 


0.09 


Arg 


R 


135.583 


1.333 


108.86 


10.28 


2.88 


2.52 


-3.44 


Asn 


N 


72.969 


-0.753 


-55.42 


6.17 


3.22 


1.45 


0.84 


Asp 


D 


68.631 


-0.220 


47.89 


3.37 


3.64 


1.13 


2.36 


Cys 


C 


65.063 


0.477 


160.13 


0.02 


0.71 


-0.97 


4.13 


Gin 


Q 


95.356 


-0.344 


134.68 


1.42 


2.18 


0.53 


-1.14 


Glu 


E 


87.730 


0.384 


53.27 


2.37 


3.08 


0.39 


-0.07 


Gly 


G 


7.206 


-0.616 


-148.03 


n.d. 


2.23 


-5.36 


0.30 


His 


H 


100.603 


0.249 


24.57 


3.38 


2.41 


1.74 


1.11 


lie 


1 


98.084 


1.529 


-104.89 


0.05 


-4.44 


-1.68 


-1.03 


Leu 


L 


102.422 


1.678 


-148.53 


2.48 


-4.19 


-1.03 


-0.98 


Lys 


K 


116.134 


1.092 


47.61 


2.78 


2.84 


1.41 


-3.14 


Met 


M 


107.039 


1.581 


6.37 


2.79 


-2.49 


-0.27 


-0.41 


Phe 


F 


121.451 


2.021 


47.67 


3.59 


-4.92 


1:30 


0.45 


Pro 


P 


77.446 


1.005 


169.73 


^ 10.49 


-1.22 


0.88 


2.23 


Ser - 


S 


44.705 


-0.733 


30.24 


0.73 


1.96 


-1.63 


0.57 


Tlir 


T 


67.162 


-0.043 


46.04 


3.71 


0.92 


-2.09 


-1.40 


Trp 


W 


153.493 


2.379 


178.69 


0.08 


-4.75 


3.65 


0.85 


Tyr 


Y 


128.937 


1.547 


49.11 


3.50 


-1.39 


2.32 


0.01 


Val 


V 


79.965 


1.117 


-106.54 


6.54 


-2.69 


-2.53 


-1.29 



Ligands interact with biological target proteins v/a various forces such as ionic 
interactions, ion-dipole interactions, dipole-dipole interactions, hydrogen bond 
5 interactions, hydrophobic interactions, 7i-stacking interactions, edge-on aromatic 
interactions, cation-7t interactions, dispersion and induction forces. Accordingly, 
physlcochemlcal descriptors reflecting these Interaction forces have successfully been 
employed in descriptors used in Quantitative Structure-Activity Relationships (QSAR), 
. Principle Component Regression (PGR) and Partial Least-Squares (PLS) analysis of 
10 drug/ligand responses. "Quantitative structure-activity relationships and experimental 
design", U. Norinder and T. Hogberg, "Textbook of drug design and discoverf, Taylor 
and Francis, London, (2002), ppl 17-1 55. 

The physicochemical descriptors can be experimentally derived and/or theoretically 
15 calculated. The descriptors can be seen to reflect hydrophobic properties, electronic 
properties, steric properties or hydrogen bonding capabilities. Some descriptors can be 
seen to reflect combinations of such properties, especially combinations of electronic 
and steric features. The present invention describes a method or methods \A^erein the 
physicochemical descriptors reflect 7TM receptor-ligand interaction features of the 
20 amino acid residues. Additionally, the physicochemical descriptors are chosen to 
reflect hydrophobic, electronic, steric, hydrogen bonding or other properties of the 
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amino acid residues. Yet further, tlie physicocfiemical descriptors may reflect 3- 
dirhensiional features of the amino acid residues. 

The physlcochemicai descriptors of the present method may be selected from 
5 descriptors used in quantitative structure-activity relationships (QSAR), Principle 
Component Regi-esslon (PCR) and Partial Least-Squares (PLS) analysis of. peptides. 

Typical hydrophobic descriptors are e.g. Partition coefficient (logP), Calculated partition 
coefficient (clog P, Prolog P, Maclog P), Distiibution coefficient (log D), Polar surface 
10 area, Nonpolar surface area, TLC retention time, HPLC retention time, and HPLC 
capacity factor (log l<). 

Typical steric parameters are e.g. Molecular weight (MW), van der Waals volume, van 
der Weals radius. Molar refractivity (MR), STERIMOL parameters (L. Bi. Bs), Total 
15 surface area, occupied volume by a residue buried in globular protein, and buikiness 
defined as tiie ratio of the side-chain volume to its length. 

Typical electronic parameters are e.g. lonisation constant (pKcooH. pKnh2), Isoelectric 
point, Net charge at pH 7, NMR chemical shitt, ^^C NMR chemical shift. Calculated 

20 interaction energies, Electronic Charge Index (ECl), Charge transfer for carbons (CT), 
Maximum electrostatic potential (Vmax). Minimum electrostatic potential (Vmih)r Maximum 
local ionization energy (Imax), Minimum local ionization energy (Imm), Molecular 
Electrostatic Potential (MEP) on Connolly Molecular Surface, Energy of highest 
occupied molecular orbital (Ehomo). Energy of lowest unoccupied molecular orbital 

25 (Elumo), Pipole moment (p), Polarizabillty (a), Most positive partial charge on a 

hydrogen atom (qH+), Most negative partial charge in the molecule (q-), and Partial 
charges on the oxygen and carbon atoms (qC, qO) of the carbonyl group. 

Thus, the physlcochemicai descriptors may be selected from molecular weight (MW), 
30 van der Waals volume, van der Waals radius, molar refractivity (MR), STERIMOL 

parameters (L, Bi, 85), Parachor (Pr), polar surface area, non-polar surface area, total 
surface area, lonisation constant (pKcooH, pKnh2), isoelectric point, net charge at pH 7, 
partition coefficient (log P), calculated partition coefficient (clog P, Prolog P, Maclog P), 
distribution coefficient (log D), TLC retention time, HPLC retention time, HPLC capacity 
35 factor log k, 'H NMR chemical shift, "C NMR chemical shift, steric and electrostatic 
3D-property MS-WHIM Indexes, calculated Interaction energies, isotropic surface area 
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(ISA), electronic charge index (ECl), charge transfer for carbons (CT), Lewis basicity 
(LB). Lewis acidity (LA), maximum electrostatic potential (Vmax). minimum electrostatic 
potential (Vmin). maximum local ionization energy (Imax). minimum local ionization energy 
(Imin). confonnational strain energy (AHstrain), molecular electrostatic potential (MEP) on 

5 Connolly molecular surface, local flexibility (Fr), flexibility index (Fb), chain flexibiiity 
(FO), occupied volume by a residue buried in globular protein, bulkiness defined as the 
ratio of the side-chain volume to its length, total energy (Etotai), heat of fomnatlon (AHf ), 
energy of highest occupied molecular orbital (Ehomo). energy of lowest unoccupied 
molecular orbital (Etuwo). dipole moment (p), polarizability (a), most positive partial 

10 charge on a hydrogen atom (qH+), most negative partial charge in the molecule (q-), 
partial charges on the oxygen and carbon atoms (qC, qO) of the carbonyl group, 
integrated molecular transform (FTm), integrated electronic transfomi (FTe), Integrated 
charge transform (FTc), normalized molecular moment (Mn), electronic moment (Me), 
charge moment (Mc), absolute electronegativity (EN), absolute hardness (HA). Such 

15 descriptors convey infonmation on ligand-binding features in a biological target protein 
such as transmembrane receptors including 7TM receptors. 

The physicochemical descriptors of amino acids or of amino acid side chains can also 
be obtained from principal component analysis (PCA) of the above-mentioned 

20 physicochemical descriptors, e.g. such as principal properties z-scales derived from 
collections of experimental data or with additional theoretical descriptors, MS-WHIM 
3D-description matrices reflecting structural and electronic features of molecules, t- 
scores from interaction energies calculated with program GRID, and other 
combinations of descriptors mentioned above. C.f. Priolo et al., J. M. Journal of 

25 Molecular Catalysis B: Enzymatic (2001) 15, 177-189, Lejon et al.. Journal of Peptide 
Science (2001) 7, 74-81 , Zaiiani et al., Journal of Chemical Infonnation and Computer 
Sciences (1999) 39, 525-533, Matter, H. Journal of Peptide Research (1998) 52, 305- 
314, Sandberg et al. Journal of Medicinal Chemistry (1998) 41, 2481-2491 , Collantes 
and Dunn, Journal of Medicinal Chemistry (1 995) 38, 2705-1 3, Norinder, Ulf . Peptides 

30 (New York. NY, United States) (1991) 12, 1223-7, Hellberg and Kern. International 
Journal of Peptides, Protein Research (1990) 36, 440-4, Skagertjerg et al. Journal of 
Chemometrics (1990) 4, 241-53, Hellberg et al. Journal of Medicinal Chemistry (1987) 
30, 1 126-1 135. In other words, In the methods according to the present invention 
wherein step v) is included, a simplified measure of the physicochemical properties of 

35 the binding site Is obtained from principal component analysis (PCA) of the 
physicochenriical descriptors. 
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Each residue type may be assigned as many physicochemicai descriptors as decided, 
providing additional detalis of cliemlcal features of the binding site of interest In the 
following, a bitmap for a given selection of binding site residues is denoted F. Other 
5 descriptors used in the references cited herein may be chosen or combinations of 
these or novel descriptors reflecting physicochemicai properties of amino acids 
relevant for ligand receptor interactions may be selected. 

The physicochemicai descriptors according to the present invention may also include 
10 dummy parameters or indicator variables, e.g. 1 and 0. Said indicator variables may 
denote the absence or the presence of aromatic side chains, hydrophobic side chains, 
negatively charged side chains, positively charged side chains, polar side chains, 
hydrogen-bond donating side chains, hydrogen-bond accepting side chains and/or 
other selected features. 

15 

The binding site classification procedure using for example binary codes is illustrated in 
Figure 3. 

In this example, a normalised string of bits, 0 or 1 representing chemical features of the 
20 binding site residues are generated. A set of five bits are assigned to each amino acid 
residue specifying the absence 0 or presence 1 of a certain chemical feature or 
characteristics. In this example the indicator variables correspond to the presence 
(TRUE = 1 ) or absence (FALSE = 0) of hydrophobic, aromatic, positively charged, 
negatively charged or polar features. For example, a tyrosine residue using such 
25 indicator variables will be represented by the bitmap fingerprint 110 01 (hydrophobic 
- aromatic - absent - absent - polar). The mapping process of physicochemicai 
descriptors into a string containing information of all selected amino acids is usually 
carried out by conventional computerized methods. Certain chemical features may be 
considered more or less important than others, and weighted accordingly in the binding 
30 site classification. The present invention therefore describes a method wherein step v) 
is included and the physicochemicai descriptors are weighted in step v). 

An embodiment of the invention uses pseudo-sequences comprising at the most 50 
amino acids obtained from at the most 12 amino acid residues per 7TM helix or 
35 extracellular loops, sequential or non-sequential, which are associated with 

physicochemicai descriptors reflecting hydrophobic, electronic, steric, and hydrogen 
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bonding properties. A specific embodiment of the invention uses pseudo-sequences 
comprising at tlie most 40 amino acids obtained from at tlie most 8 amino acid 
residues per 7TM helix or extracellular loops, sequential or hon^sequential, which are 
associated with physicochemical descriptors reflecting hydrophobic, electronic, steric, 
5 and hydrogen bonding properties. A prefen-ed embodinient of the invention uses 

pseudo-sequences comprising at the most 30 amino acids obtained from at the most 6 
amino acid residues per 7TM helix or extracellular loops, sequential or non-sequential, 
which are associated with theoretically derived physicochemical descriptors reflecting 
hydrophobic, electronic, steric, and hydrogen bonding properties. 
■10- " ■ 

Similarity Scores 

It is desirable to quantify how similar a given receptor binding site or subsite is to other 
receptor binding sites and or subsites. Here we apply a number of different similarity 
measures to ranl< binding sites and their corresponding bitmaps. The measures are 

15 chosen due to their capabilities to handle large data sets and to iteratively, if heeded, 
allow the investigation of relationships between combinations of different subsites or 
binding sites. The measures could handle different types of descriptors described 
herein and may be based upon a pattern recognition method, a Principal Component 
Analysis (PCA) reducing the number of descriptors to a few principal components, a 

20 Tanimoto Similarity Measure, a Tversky Similarity Measure or a Euclidian Distance 
Measure as described in Press, W.H; Flannery, B.P.; Teukolsky, S.A.; Vettering, W.T. 
Numerical recipes: The art: of scientific ooniputing; Cambridge University Press: 1986. 
The present invention relates to methods, wherein the generation of a similarity score 
in step v) is based upon a pattern recognition method. In an embodiment of the 

25 methods according to the present invention, the generation of the similarity score 

involves a Principal Component Analysis (PCA) reducing the number of descriptors to 
a few principal components. 

The similarity measure applied are the Tanimoto Coefficient TO, Tversky similarity 7S, 
30 and an Euclidian distance d(F1,F2) defined below. 

Tanimoto Similarity Measure: 

The Tanimoto coefficient between two bitmaps F1 and F2 is defined as 



35 



TC-BC/(B1*B2-BC) 
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Here B1 and B2 are the numbers of 1's \r\ F1 and F2 respectively and BC is the 
number of 1's In common between F1 and FZ TC is a value between 0 and 1. If TC = 1 
F1 and F2 are identical. If TC = 0, F1 and F2 have no TRUE = 1 occasions in common. 
In a further embodiment of the methods according to the present invention the 
5 generation of the similarity score in step v) (above) is based upon a Tanimoto Similarity 
Measure: tC = BG / (B1 + B2 - BC). 

Tversky Similarity Measure: 

The Tversky coefficient between two bitmaps F1 and F2 Is defined as 

10 ■•. 

^ TSrBC/(a*B1Unique'^fi*B2Unique'¥BC) ^ 

Here 81 Unique and B2Unique are the number of unique 1 's in F1 and F2 respectively. 
a and fi are constants used to weight prototype and variant features. When a = >9 = 1 , 
15 this measure produces a symmetrical similarity, metric identical to TC. In a still further 
embodiment of the methods according to the present invention, the generation of the 
similarity score In step v) is based upon a Tversky Similarity Measure: TC = BC / (a * 
81 Unique + p * B2Unique + BC), wherein a are prototype features and p variant 
features. 

20 

Euclidian Distance Measure 

The Euclidian distance represents the geometric distance between the bitmaps F1 and 
F2 

25 d(F1,F2)-sqrt(F1'F2)^ 

where F1 and F2 are vectors in a N dimensional space. In the situation where 
physicochemical descriptors are used, F1 and F2 are no longer bit strings but string 
containing physicochemical descriptors representing the binding site residues of 
30 interest. In yet another embodiment of the methods according to the present Invention, 
the generation of the similarity score in step v) is based upon Euclidian Distance 
Measure: d(F1 .F2) = sqrt (F1 ~ F2)l 

Ranking of 7TI\/I receptors 
35 A ranking of the 7TM receptors based on the physicochemical properties of their 

binding sites may be obtained, which gives a good indication of the similarity between 
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them. In certain cases thjs is important, and in such cases,, step vl) (above) is included. 



The ranking with respect to the physlcochemlca! properties assigned to the aligned 
pseudo-sequence Is based upon similarity scores obtained according to the procedures 
5 described aboye, or from distances between coordinates In a Principle Component 
(PC) n-dimenslonal space, but it may also be based upon a 2- or 3-dimenslonal 
graphical representation. In the latter case, visualisation of the relationship and 
similarities between 7TM receptors Is simplified. 

10 Examples of receptors where sequence alignment and the physicogenomics approach 
give comparable relationships between 7TM receptors are typically illustrated by 
certain receptors with subclasses such as the neurokinin NK1 to UK4 receptors and 
muscarinic M1 to M5. 

15 Thus, a model using theoretical descriptors will rank the receptors closest to the 
muscarinic MZ as M5 (3.4), M2 (3.7). Ml (3.9) and M4 (4.3), which is in accordance 
with findings that muscarinic antagonists are in principle fairly subtype-unselective. 
Likewise, the same model ranks the receptors closest to the neurokinin NK1 as NK3 

(3.3) , NK4 (3.3) and NK2 (3,6), followed by Adenosine A3R (3.6). 

20 

in contrast, the same model applied to histamine H2 will rank the closest receptors as 
adrenergic b1 (3.6) and b3(3.7), whereas the closest histamine receptor is more 
remote, I.e. histamine HI (4.3) and the remaining ones are even further away but very 
close to each other. I.e. 1-13 (5.5) and H4 (5.5). 

25 

The same model applied to GRP44 (CRTH2) discussed previously will rank the closest 
receptors as angiotensin AT2 (2.9), chemoklne-like receptor 1 (3.1 ), bradykinin B2 

(3.4) , Chemokine Receptor type 10 (3.6), Chemokine Receptor 1 (3.6), and angiotensin 
ATI (3.7), which where not identified in the phylogenetic tree as close neighbours. 

30 

In one embodiment, the amino acid residues, up to six per helix, are selected from TM- 
III, TM-IV, TM-V, TM-VI and TM-VIII to form the following pseudo-sequences, which 
are used in the alignment. The following rank order of the similarity of the receptors can 
be obtained by implying the given set of amino adds associated with theoretically 
35 derived physicochemical descriptors reflecting hydrophobic, electronic, steric, and 
hydrogen bonding properties: 
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Receptor. Pseudosequence Ranking 

GP44_HUiyiAN HSFFMFNTYAKFAWYHSEALTA 1 

AG22_HUMAN FGLTMFSSTAKNGWFHTDALIG 2 

CML1_HUMAN SNLIMFLSSTRFGWYHNELLTA 3 

BRB2_HUMAN VNISLYLSMNLNGWFQTDTTSA 4 

CKRA_HUMAN ISYSFHLAAAQVGQYSLDTLSA 5 

C3X1_HUMAN TTFFFFVAQNTNGWYNIETLEA 6 

AG2R_HUMAN ASVSLYASAGKNGWHQTDVMIA 7 



The use of lead structures might be based on /n silico searches of specific scaffolds 
identified, on a pharmacophore model derived from these compounds, on in silico 

5 . searches of such pharmacophore model, on design and synthesis of chemical libraries 
encompassing specific scaffolds identified, on a phamnacophore model derived from 
these compounds to design and/or construct chemical libraries containing novel 
chemical features compatible with the pharmacophores, on a pharmacophore models 
derived from these compounds to specifically design and synthesise novel ligands or 

10 on other common technologies used in drug design. The present invention also relates 
to the use of a pharrnacophore as described herein for in silico screening, for 
construction of a library or for design of a llgand. 

Having compared and ranked the 7TM receptors to each other by a suitable 
15 computerised mathematical model based on the relevant binding site and associated 
physicochemlcal descriptors of amino acid residues occupying the binding site, one 
can also identify 7TM receptors with binding properties similar to a given receptor and 
use that infonnation to facilitate the drug design process. 

20 For example the methods of the present invention allow identification of receptors, 
which are likely to cause a selectivity problem during drug development of a drug 
interacting with a given receptor. These potentially interfering receptors could be 
subject to directed counter-screens, reducing the need to screen compounds very 
broadly on a large number of receptors in the dmg discovery and development 

25 process. The present invention relates to the use of the methods described herein to 
Identify receptors, which are likely to cause a selectivity problem during drug 
development of a drug Interacting with a given receptor. 
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Analogously to the utilization of the method to Identify similarities in binding sites, the 
same principle can be applied to identify differences in subsites of binding sites 
between 7TM receptors as means to improve receptor selectivity of a drug towards a 
given TTM receptor. The information regarding where significant differences exist in 
5 subsites between related receptors can be used in the design of ligands with improved 
receptor selectivity of a drug towards a given TTM receptor. 

Legends to figures 

Figure 1 illustrates the conventional phylogenetic analysis of the GPR44 {CRTH2) 
10 receptor, 

Figure 2.1 shows a schematic depiction of the secondary structure of a rhodopsin-like 
7TM receptor with one or two conserved, key residues highlighted in each 
transmembrane segment: Asnl:18; Aspll:10; Cyslll:01 and Arglll;26; TrplV:10; ProV:16; 
15 ProVI:15; ProV!l:17, 



Figure 2.2 shows how the transmembrane segments are generically numbered based 
on of the key residues present in the family B class of receptors, 

20 Figure 3 shows the binary 5-digit codes used to indicate absence or presence 6f 
physicochemlcal descriptor. 

The following example is included for illustrative purposes and is not intended to limit 
- the invention in any way. 



25 
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Example 

A moder using theoretical size and electronic descriptors ranked the angiotensin AT2 
and angiotensin AT1 receptors to be closely related to GRP44 (CRTH2) with respect to 
ligand-blnding features in the binding site. Accordingly, among identified ligands, the 
i<nown ATI antagonist candesartan was found to Inhibit [^H]PGD2 of GRP44 with an 
ICeoOf 2.1 ^iM. 




CANDESARTAN 



10 This example - describing the alignment procedure followed by amino acid selection 
and the other steps outlined above leading to the identification of ATI as a GRP44- 
related receptor that serve as a source of discovering lead molecules for the new target 
GRP44 that Can be used for In silico screening or design of focused chemical libraries - 
can arialqgously be applied to other receptors 0 

15 



PCT/PK2004/000148 

33 

1. A pseudo-sequence method for comparing a first 7TM receptor with one or more 
further 7TM receptors with respect to the physicochemical properties of selected amino 

5 acid residues of their binding sites, the method comprising the steps of: 

i) optionally, aligning part of or all of the amino acid sequence of the first 7TM receptor 
with part of or all of the amino acid sequence of the one or more further 7TM receptors, 

10 ii) selecting, in a sequential or non-sequential order, at the most 1 2 amino acid 

residues per helix and/pr extracellular loops, which are involved in one or more binding 
sites of each 7TM receptor, 

iii) forming a pseudo-sequence comprising at the most 50 amino acid residues frorn the 
15 selected sequential or non-sequential amino acid residues, 

Iv) for each 7TM receptor assigning one or more physicochemical descriptors to the 
amino acid residues of the selected amino acid pseudo-sequence involved in one or 
more binding sites, 

20 

v) optionally, for each 7TM receptor mathematically manipulating the physicochemical 
descriptors of step iv) to obtain a simplified measure of the physjcochemical properties 
of the binding site, 

25 vi) for each 7TM receptor generating a similarity score as defined herein by comparing 
the physicochemical descriptor or, if relevant, the simplified measure for the first 7TM 
receptor with the physicochemical descriptors or, if relevant, the simplified measures 
for the one or further 7TM receptors, 

30 vii) optionally, ranking the 7TM receptors with respect to the physicochemical 
properties of their binding sites according to the similarity scores obtained in step vi). 

2. A method according to claim 1, wherein the comparison is made without using data 
related to binding affinity of a ligand to a 7TM receptor. 

35 

3. A method according to claim 1 for classifying 7TM receptors according to the 
physicochemical properties of their binding sites. 
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4. A method according to claim 3, wherein the classification is made without using data 
related to binding affinity of a llgand to a 7TM receptor. 

5 , 5. A method according to any of the preceding claims, \wherein step 11) as defined In 
claim 1 comprises selecting, In a sequential or non-sequential order, at the most 11 
such as, e.g., at the most 10, at the most 9, at this most 8, at the most 7 or at the most 
6 amino acid residues per helix and/or extracellular loops, which are Involved in one or 
more binding sites of each 7TM receptor. 

10 

6. A method according to any of the preceding claims, wherein step ill) as defined in 
claim 1 comprises forming a pseudo-sequence comprising at the most 50 such as, e.g., 
at the most 45, at the most 40, at the most 35 or at the most 30 amino acid residues 
from the selected sequential or non-sequential amino add residues. 

15 ■ . . . 

7. A drug discovery method for identifying ligands, which b|nd to a first 7TM receptor 
and potentially bind to one or more further 7TM receptors, the method comprising the 
steps of I) to vll) as defined In claim 1 , 5 or 6 and the further steps of 

20 viii) selecting from one to about 100 further 7TM receptors which have the closest 
similarity scores to the first 7TM receptor, 

ix) identifying ligands which potentially bind to those further 7TM receptors selected In 
step vll) by selecting ligands that bind to the first 7tM receptor. 

25. 

8. A drug discovery method for identifying ligands which bind to a first 7TM receptor 
and to one or more further 7TM receptors, the method comprising the steps of I) to vll) 
as defined in claim 1 , 5 or 6 and the further steps of: 

30 viii) selecting from one to about 100 further 7TM receptors which have the closest 
similarity scores to the first 7TM receptor, 

ix) screening ligands that bind to the first 7TM receptor against the selected 7TM 
receptors of step viii). 

35 

9. A drug discovery method for identifying a potential lead compound for a first TTIVI 
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receptor, the method cx)mprising the steps of I) to vli) as defined in claim 1, 5 or 6 and 
the further steps of 

viil) selecting from one to about 100 further 7TM receptors which have the closest 
5 similarity scores to the first 7TM receptor, 

ix) Identifying ligands that bind to said one or more further TTIVl receptors to construct a 
library including a potential lead compound for the first 7TM receptor. 

10 10. A drug discovery method for identifying a lead compound for a first 7TM receptor, 
the method comprising the steps of i) to vii) as defined in claim 1 , 5 or 6 and the further 
steps of 

viii) selecting from one to about 100 further 7TM receptors which have the closest 
15 similarity scores to the first 7TI\/I receptor, 

ix) identifying ligands that bind to said one or more further 7TM receptors to construct a' 
library, and 

20 x) screening said library against the first 7TM receptor to identify a lead compound for 
the first 7TM receptor. 

1 1 . A drug discovery method for constructing a pharmacophore model for a first 7TM 
- receptor, the method comprising the steps of i) to vii) as defined in claim 1 , 5 or 6 and 

25 the further steps of 

viii) selecting from one to about 100 further 7TM receptors which have the closest 
similarity scores to the first 7TM receptor, 

30 ix) Identifying ligands that bind to said one or more further 7TM receptors to construct a 
pharmacophore model. 

12. A drug discovery method according to claim 10, wherein the first 7TM receptor is 
one for which no ligands have been identified. 

35 
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1 3. A drug discovery method according to claim 1 0 or 1 1 , wherein the first 7TM 
receptor is an orphan receptor. 

14. A method according to any of claims 7-12, wherein from one to 50 further 7TM 
5 receptors is/are selected in step viii). 

1 5. A method according to any of claims 7-12, wherein from one to 25 further 7TM 
receptors is/are selected in step viii). 

10 16. A method according to any of claims 7-12, wherein from one to 15 further 7TM 
receptors is/are selected in step viii). . ^ 

17. A method according to any of the preceding claims, wherein the method is 
executed by a computer under the control of a program and the computer includes a 

15 memory for storing said program. 

1 8. A method according to any of the preceding claims, wherein step i) is included and 
the alignment is based on a model developed for 7TM receptors. 

20 19. A method according to claim 18, wherein the 7TM receptors are Class A, Class B, 
Class C or taste receptors. 

20. A method according to any of the preceding claims, wherein step i) is included and 
the alignment is made with respect to transmembrane positioning of a-helices of 7TM 

25 receptors. 

21 . A method according to any of the preceding claims, wherein the binding site 
includes amino acid residues located in one or more extracellular loops of the 7TM 
receptors. 



30 



22. A method according to any of the preceding claims, wherein the binding site 
includes amino acid residues located In one or more subsites of the binding site and in 
one or more extracellular loops of the 7TM receptors. 



35 



23. A method according to any of the preceding claims, wherein the physicbchemical 
descriptors reflect 7TM receptor-ligand interaction features of the amino acid residues. 
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24; A method according to any of the preceding claims, wherein the physicocheniical 
. descriptors are chosen to reflect hydrophobic, electronic, steric, hydrogen bonding or 
other properties of the amino acid residues. 

5 ■ ■ 

25i A method acdording to any of the preceding claims, wherein the physicochemlcal 
descriptors reflect 3-dimensional features of the amino acid residues. 

26. A method according to any of the preceding claims, wherein the physicochemical 
10 descriptors are selected from descriptors used in quantitative structure-activity 

relationships (QSAR), Principle Component Regression (PGR) and Partial Least- 
Squares (PLS) analysis of peptides. 

27. A method according to any of claims 23-26, wherein the physicochemical 

1 5 descriptors are selected from molecular weight (MW), van der Waals volume, van der 
Waals radius, molar refractivity (MR), STERIMOL parameters (L, Bi , Bg), Parachor (Pr), 
polar surface area, non-polar surface area, total surface area, ionisation constant 
(pKcooH, pKnh2). Isoelectric point, net charge at pH 7, partition coefficient (log P), 
calculated partition coefficient (clog P, Prolog P, Madog P), distribution coefficient (log 

20 D), TLC retention time, HPLC retention time, HPLC capacity factor log k, NMR 
chemical shift, "C NMR chemical shift, steric and electrostatic 3D-property MS-WHIM 
indexes, calculated interaction energies, isotropic surface area (ISA), electronic charge 
index (ECl), charge transfer for carbons (CT), Lewis basicity (LB), Lewis acidity (LA), 
rriaximum ejectrostatic potential (Vmax). minimum electrostatic potentlal (V^in), maximum 

25 local ionization energy (Lax), minimum local Ionization energy (Un), conformational 
strain energy (AHstrain), molecular electrostatic potential (MEP) on Connolly molecular 
surface, local flexibility (Fr), flexibility index (Fb), chain flexibility (FO), occupied volume 
by a residue buried in globular protein, bulkiness defined as the ratio of the side-chain 
volume to its length, total energy (Eiotai), heat of formation (AHf ), energy of highest 

30 occupied molecular orbital (Ehomo)i energy of lowest unoccupied molecular orbital 
(Elumo). dipole moment (p), polarizability {a), most positive partial charge on a 
hydrogen atom (qH+), most negative partial charge in the molecule (q-), partial 
charges on the oxygen and carbon atoms (qC, qO) of the carbonyl group, integrated 
molecular transform (FTm), integrated electronic transfonn (FJe), Integrated charge 

35 transfonrn (FTc), nomialized molecular moment (Mn), electronic mornent (Me), charge 
moment (Mc), absolute electronegativity (EN), absolute hardness (HA). 
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28. A method according to any of claims 17-25, wherein the physicochemical 
descriptors include Indicator variables such as, e.g., 1 and 0. 

5 29. A method according to claim 28, wherein the indicator variables denote absence or 
presence of aromatic side chains, hydrophobic side chains, negatively charged side 
chains, positively charged side chains, polar side chains, hydrogen-bond donating side 
chains, hydrogen-bond accepting side chains and/or other selected features. 



10 30. A method according to any of the preceding claims including step v), wherein the 
physicochemical descriptors are weighted in step v). 

31 . A method according to any of the preceding claims including step v), wherein a 
simplified measure of the physicochemical properties of the binding site is obtained 

15 from principal component analysis (PCA) of the physicochemical descriptors. 

32. A method according to any of the preceding claims, wherein the generation of a 
similarity score in step vi) is based upon a pattern recognition method. 

20 33. A method according to any of the preceding claims, wherein the generation of the 
similarity scor0 involves a Principal Component Analysis (PCA) reducing the number of 
descriptors to a few principal components. 

^34^ A niethod according to any^of the preceding claims, wherein the generation of the 

25 similarity score in step vi) is based upon Euclidian Distance Measure: d(F1 ,F2) = sqrt 
(F1-F2)l 

35. A method according to claim 28 or 29, wherein the generation of the similarity score 
in step v) is based upon a Tanimoto Similarity Measure: TC = BC / (B1 + 82 - BC). 

30 

36. A method according to claim 28 or 29, wherein the generation of the similarity score 
in step V) is based upon a Tversky Similarity Measure: TC = BC / (a * B1 Unique + p * 
B2Unique + BC), wherein a are prototype features and p variant features. 
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37. A method according to any of the preceding claims, wherein the step vii) is 
included. 
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38. A method according to any of the preceding claims, wherein the similarity score or, 
If relevant, the ranking is based upon a 2- or 3-dimensional graphical represisntation. 

5 39. Use of a phanmacophore according to claim 1 1 for in silico screening, 

40. Use of a phamiacophore according to claim 1 1 for constmction of a library. 

41. Use of a pharmacophore according to claim 1 1 for design of a ligand. 

10 

42. Use of a method according to any of claims 1-38 to identify receptors, which are 
likely to cause a selectivity problem during drug development of a drug interacting with 
a given receptor. 

15 43. Use of a method according to any of claims 1-38 to identify differences in subsites 
of binding sites between 7TM receptors as means to improve receptor selectivity of a 
drug towards a given 7TM receptor. 
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